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ABSTRACT 


In  this  research  a  general  purpose  digitally  controlled 
analog  filter  is  presented.  The  novel  design  is  a  cascade  of 
second  order  sections  that  are  individually  programmed  to 
achieve  any  filtering  topologies.  Two-binary  words  are  used 
to  control  the  pole  frequency  co  p  and  selectivity  Qp  of  each 
section  independently.  Each  second-order  section 
is  a  Generalized-Immittance  Converter  (GIC)  biguads  which  are 
known  for  their  high  stability  and  low  active  and  passive 
sensitivity.  CMOS  switches  are  used  to  electronically 
relocate  the  minimum  number  of  passive  elements  to  achieve 
function  programmability.  Switches  are  also  used  to  select 
the  number  of  cascaded  sections  to  realize  higher  order 
transfer  functions. 
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the  filter  function  with  linear  components  yet  retain  the 
flexibility  of  varying  the  filter  parameters  digitally. 

Recently,  the  advantages  of  combining  linear  components 
(operational  amplifiers,  resistors  and  capacitors)  and  non 
linear  elements  (switches)  have  been  demonstrated  using 
switched  capacitor  techniques  [36-38].  In  this  research,  we 
are  presenting  the  results  of  realizing  a  continuous  active 
device  using  linear  elements  and  switches  controlled  by 
digital  signals  to  achieve  a  fully  programmable  filter  [32]. 
Several  programming  features  of  the  proposed  filter  are 
reported.  The  first  feature  is  the  ability  of  the  networ)c  to 
realize  the  most  common  filtering  functions  (function 
programmability)  namely:  Low  Pass  (LP),  Band  Pass  (BP),  High 
Pass  (HP),  All  Pass  (AP)  and  Notch  (N)  functions,  using  the 
minimal  set  of  elements.  The  second  feature  is  the  ability 
of  the  network  to  program  (independently)  the  key  parameters 
of  the  filtering  function  chosen  (parameter  programmability) 
namely:  the  pole  resonant  frequency  (u>p)  and  selectivity 
(Qp).  Finally,  the  ability  to  program  the  network  to  cascade 
several  sections  to  achieve  higher  order  filer.  All  of  the 
above  programmability  features  are  performed  independently  to 
realize  a  universal  filtering  network.  In  order  to 
demonstrate  the  idea  of  this  research,  it  is  necessary  to 
introduce  some  theoretical  back  ground. 
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B.  LINEAR  SYSTEMS 


K 


The  box  at  Fig.  (1.1)  illustrates  the  concept  of  the 
linear  system.  In  the  time  domain,  the  system  is 
characterized  by  its  impulse  response  y(t),  which  is  the 
output  signal  y(t)  produced  in  response  to  an  impulse,  for 
an  arbitrary  input  signal  u(t),  the  output  signal  y(t)  is 
given  by  the  well  known  convolution  integral. 


% 


uLl-T')  -kLi) 
d 


Ult) 

1  ulO 

_ Fig.  1.1  -  Black  Box  Concept  of  Linear  System 

The  system  transfer  function  G(s)  is  the  Laplace 
transform  of  g(t),  thus: 


Ql’3\  -  ClI  cxp  C-5  0 
C.  S  =  iT  \  uo 


(1.2) 


* 


Under  the  laplace  transformation  1.1)  becomes: 


so  that  the  transfer  function  G(s)  is  the  ratio  of  the  output 
variable  to  the  input  variable. 


(1.4) 


The  most  general  types  of  linear  systems  are  consisting 
of  a  finite  number  of  lumped,  linear  and  time-invariant 
elements.  The  system  is  characterized  by  an  Nth  order 
ordinary  linear  differential  equation  which  results,  in  most 
general  cases,  a  transfer  function  G(s)  which  is  a  real 
rational  polynomial  function  of  the  complex  variable  S.  Thus 
we  can  write  Q(s)  in  general  as 


'  Lh  i\  U-.  s"".  V 1. 


(1.5) 


where  Pk  and  Ik  are  real  numbers  so  that  G(s)  is  real  for 

real  s,  and  the  roots  of  the  polynomial  P(s)  and  E(s)  must  be 
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where  C  is  a  constant  extracted  from  G(s)  such  that  E(s)  and 
P(s)  become  monic  polynomials  (leading  coefficients  equal  to 
unity) ,  Szk  are  the  transmission  zeros  and  Snk  are  the 
natural  modes  at  the  system. 

C.  FILTERS  AS  A  SPECIAL  CLASS  OF  LINEAR  SYSTEMS 

Linear  systems  can  be  distinguished  into  "SPECTRAL 
SHAPING  NETWORKS"  and  "FILTERS."  The  role  of  filters  is  one 
of  selecting  signals  while  the  role  of  spectral  shaping 
networks  is  that  of  modifying  the  input  signal  spectrum  in  an 
arbitrary,  but  predescribed  manner.  Specifically,  we  desire 
that  a  filter  should  do  as  little  as  possible  shaping  on 
signals  in  its  passband;  any  shaping  it  is  considered  a 
distortion  of  the  signal.  On  the  other  hand,  networks  which 
perform  pulse  forming  fall  within  the  spectral  shaping 
category. 

D.  ACTIVE  FILTER  FUNDAMENTALS 

The  distinction  between  passive  and  active  filters  is 
that  the  first  do  not  require  a  power  source  to  perform  their 
function  while  the  second  do.  The  motivation  behind  active 


RC  filters  lies  in  the  desire  to  have  inductorless  filter 


realizations.  It  is  well  known  that  of  the  three  passive 
R,C&L  elements,  the  inductor  is  the  most  nonideal  one.  This 
is  especially  true  at  low  frequencies,  where  inductors  become 
quite  bully  and  have  increased  losses  or  equivalently  lower 
Q-f actors. 

E.  GENERALIZED  IMMITTANCE  CONVERTOR 

One  of  the  methods  of  active  RC  filters  design  consists 
of  simulating  the  inductances  in  the  LC  ladder  by  active  RC 
networks.  This  simulation  can  be  based  on  the  principle  that 
we  want  to  find  a  one  port  network  having  an  input  impedance. 

Zll  =  s*L 

Various  active  elements  as  well  as  synthesis  procedures 
employing  them  have  been  proposed  [1-9].  A  partial  list 
includes: 


(1) 

Negative  Impedance  Converter 

(NIC) 

(2) 

Negative  Impedance  Inverter 

(NIV) 

(3) 

Postive  Impedance  Converter 

(PIC) 

(4) 

Gyrator 

(5) 

Generalized  Impedance  Converter  (GIC) 

(6) 

Curent  Conveyor  ( CC ) , and 

• 

(7) 

Operational  Amplifier  (OA) 

Although  the  introduction  of  these  elements  has 
stimulated  research  in  the  area  of  active  network  theory, 
very  few  elements  have  made  their  way  to  large  scale 


production  to  become  available  as  off-the-shelf  items.  The 
reason  for  this  is  mostly  an  economic  one.  For  a  device  to 
become  available  at  low  cost,  it  has  to  be  used  in 
substantially  large  quantities.  It  follows  that  such  a 
device  has  to  be  versatile  enough  to  be  of  use  in  a  number  of 
applications,  of  which  active  filter  design  is  only  one.  The 
analog  circuit  design  area  has  found  these  attributes  in  the 
IC  operational  amplifiers  (Op.  Amp.) 

The  IC  Op.  Amp.  is  currently  the  most  popular  linear 
active  element.  It  is  available  from  a  large  number  of 
manufacturers,  at  reasonable  cost  and  with  good  performance 
characteristics.  Furthermore,  elements  engineers  have  become 
accustomed  to  the  use  of  Op.  Amp.  It  is  therefore,  only 
natural  that  the  Op.  Amp.  is  becoming  the  most  popular  active 
element  in  the  design  of  active  RC  filters,  and  can  be  found 
in  NIC's,  PIC's,  GIC's,  and  other  circuit  realizations. 


THEORETICAL  ANALYSIS 


A.  BIQUADRATIC  TRANSFER  FUNCTIONS 

The  filter  as  a  special  class  of  linear  system  has  a 
transfer  function  expressed  in  a  polynomial  quotient  form 
given  as 


IPkS^ 

K-O _ 

IC 


(2.1) 


where  Pk  and  gk  are  real  numbers  to  that  T(s)  is  real  for  s, 
and  the  roots  of  the  polynomials  P(s)  and  E(s)  must  either  be 
real  or  occur  in  conjugate  pairs.  Also,  in  general,  the 
degree  of  the  numerator  (deg.[P(s)]  =  M)  is  less  than  or 
equal  to  the  degree  of  denominator  (deg[E(s)]  =  N)  and  the 
roots  of  E(s)  are  in  the  open-half  S-plane.  The  E(s)  is 
known  as  the  characteristic  polynomial  or  natural  mode 
polynomial  of  the  linear  system,  and  the  degree  of  E(s),  that 
is  N,  is  the  order  or  degree  of  the  system. 

A  general  second-order  transfer  function  or  "biquad" 
function  may  be  written  as 


.«^Y*rv 


where  P(s)  is  the  loss-pole,  or  more  appropriately  here,  the 
transmission- zero  polynomial,  and  E(s)  is  the  natural  pole 
polynomial  mode  as  discussed  above.  It  is  a  usual  practice 
to  express  the  denominator  in  terms  of  _p  and  Qp,  where  _p  is 
the  natural-mode  or  resonance  frequency  and  Qp  is  the 
natural-mode  or  quality  factor.  Thus  (2.2)  becomes 


1  S  P<\.S  I-  to 


(2.3) 


The  numerator  coefficients  determine  the  location  of  the 
transmission  zeros  and  hence,  the  type  of  filter  function  the 
biquad  provides.  Special  cases  of  interest  are: 


1.  Low  Pass  (LPl 


For  which  Pl=:P2=0,  thus  two  transmission  zeros  are  at 


infinity 


5.  All  Pass  (AP) 

For  which  the  pair  of  zeros  are  at  the  mirror  image 
location  of  the  pair  of  poles,  that  is 


B.  SENSITIVITY  FUNCTIONS 

A  concern  about  the  design  of  a  filter  is  how  close  the 
resulting  response  will  be  to  the  ideal  or  desired  function. 
The  reason  for  response  deviation  from  the  ideal  is  the 
finite  tolerances  of  the  RC  elements,  as  well  as  the  nonideal 
performance  of  the  active  elements.  In  the  latter  case,  not 
only  the  gain  changes  or  tolerances  have  to  be  considered  but 
the  effect  of  the  "limited  amplifier  bandwidth"  on  the  filter 
response  must  also  be  evaluated.  Although  effects  of  initial 
component  tolerances  may  be  "trimmed  out"  during  the  initial 
filter  alignments  or  tuning  process,  a  sensitive  design  will 
deviate  from  the  required  specifications  as  time  process,  due 
to  component  variations  with  temperature,  aging,  humidity, 
etc.  Note  also  that  a  sensitive  design  might  be  extremely 
difficult  to  tune  in  the  first  place,  or  the  initial 
adjustment  will  be  quite  uneconomical. 

The  answer  to  tolerance  question  can  be  obtained  through 
sensitivity  studies.  Considerable  emphasis  has  been  placed, 
in  the  active  filter  literature,  on  the  study  of  sensitivity 


functions  and  relations,  some  of  the  most  useful  and  widely 
accepted  sensitivity  functions  are  the: 

1.  Magnitude  Function  Deviation 

Assure  a  filter  designed  to  meet  a  certain  magnitude 
characteristics  T(s)  or  T(j«^).  One  is  concerned  with  the 
deviation  in  |T(jO)j,  that  is/i.|T{  jvO  i  both  in  passband  and  in 
stopband.  Usually  it  is  desirable  to  express  the  expected 

-deAtiat ton. -  in. . dE  _ The de.viat ion  D  ( lO )  dB  in  the  magni tude 

function  may  be  evaluated  as  follows.  Let  the  function  |T(j 
change  to  [jT(j.x^)|  ^  |T(j  )|],  then 


r  QO  locj 

ICjuo> 


C^h>) 


(2.9) 


or 


Dluj:)  -  Q  bz  U 


i  + 


A  I  ’  (JVjO) 
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,(db) 


(2.10) 


and  for  small  variability  (2.10)  can  be  approximated  as 


Thus,  the  deviation  in  the  magnitude  response  in  nepers 
is  equal  to  the  per  unit  variability  in  the  magnitude  of  the 
transfer  function.  The  problem  now  reduces  to  that  of 
evaluating  the  per  unit  change  in  |T(jvP)|.  This  is  not  an 
easy  problem  since  T(joJ)  is  a  function  of  many  elements  with 
different  tolerances  and  tolerance  statistics.  Furthermore, 
the  per  unit  change  is  function  of  frequency. 

2.  Classical  Sensitivity 

Lets  recall  the  definition  of  the  classical 
sensitivity,  where  y  is  a  variable  of  interest,  usually  a 
function  of  many  parameters  of  which  x  is  one,  then 


(2.12) 


Note  that  from  the  above  definition,  is  the  limit  to  as 
Dx  — >0.  thus,  for  small  variations. 


.  i 

’  ■  j 


(2.13) 


The  usefulness  of  the  classical  sensitivity  function 
is  evident  from  (2.13).  The  per  unit  or  percentage  change  in 


Y,  due  to  a  given  per  unit  or  percentage  change  in  x,  can  be 
easily  obtained  by  multiplication  with  S,  i.e.. 


(2.14) 


3 .  Gain  Sensitivity  Product 

An  important  consideration  in  the  evaluation  of  the 
sensitivity  of  a  filter  parameter  as  considered  in  [2]  with 
respect  to  the  closed  loop  gain  is  the  tolerance  on  the 
closed  loop  due  to  the  open  loop  gain  variability. 

Thus,  the  gain-sensitivity-product,  G.S^  is  defined 
as: 


^  5, 


(2.15) 


where  )c  is  the  closed  loop  gain.  We  can  extend  (2.15)  to  the 
open  loop  gain  as: 


(2.16) 


and  also  we  can  note  that: 


(2.17) 


and  thus  then  ultimate  good  is  the  variability  or  tolerance 
rather  than  the  sensitivity,  the  gain-sensitivity  product  is 
a  better  index  for  comparing  different  designs. 

4  Determining  the  Variability  of  the  Transfer  Function 


Amplitude 

Assure  that  the  active  filter  has  i  resistors, 
rri  capacitors  and  n  amplifiers.  Let  the  amplifier  k  have  an 
open  loop  gain  Aok  and,  possibly,  a  closed  loop  gain  Kk.  The 
variability  of  the  magnitude  function  as  given  by  the 


Reference  [10]  is: 


1 

1 


•.V'.  (“-c) 


I  Ac^: 


(2.18) 


Note  that  each  of  the  sensitivity  functions  is  (2.13) 
is  a  function  of  frequency.  In  a  high  order  filter 
realization,  the  different  sensitivity  functions  might  be 
difficult  to  evaluate. 


C.  PROPOSED  GIG  FILTER  ANALYSIS 


In  order  to  obtain  the  transfer  functions  of  the  proposed 
programmable  filter  [13]  shown  at  Fig.  (2.1),  nodal  analysis 
was  used  as  follows: 

(1)  The  circuit  of  Fig.  (2.1)  was  replaced  by  the  one  of 
Fig.  (2.2)  in  which  the  two  operational 

amplifiers  were  replaced  by  the  two  equivalent  dependent 
voltage  sources  as  can  be  seen,  every  element  (node, 
admittance,  voltage  source,  etc.)  was  labeled  and  every 
element  was  associated  to  a  current  direction  and  a  voltage 
polarity. 

(2)  The  kirchoff  current  low  was  written  for  every  node 
except : 

(a)  The  reference, 

(b)  any  node  connected  to  the  reference  by  a  voltage 
source . 

.  Node  K.C.L. 

1  i7=-i3  (2.19) 

2  i2=i5+i6  (2.20) 

3 

4 


5 


i4=i7+i8 


(2.21) 


( 3 )  Every  admittance  current  was  expressed  in  terms  of 
nodal  voltages: 


il  =  yi[V3-Vl] 

(2.22) 

i2  =  Y2h/3-V2] 

(2.23) 

i3  =  Y3['v/4-Vl] 

(2.24) 

i4  =  Y4Et/4-V^] 

(2.25) 

i5  =  Y5[t/2-V6] 

(2.26) 

i6  =  Y6[V2] 

(2.27) 

i7  =[Y7[V5-V6] 

(2.28) 

i8  =Y8[ij5] 

(2.29) 

( 4 )  The  source  dependencies  were  listed  expressed  in  terms 
of  nodal  voltages. 


and  then  (substitute)  where  necessary. 

(5)  A  matrix  equation  having  the  unknown  voltages  V3,  V4, 
and  V5  related  with  the  desired  transfer  functions  Tl,  T2, 
T3,  of  the  filter  was  obtained  as: 
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(2.32) 
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and  by  substituting  ,  Tar"}!! 

'\)iit  'Unl, 


and  T3  - 


-06 


"U  i/u 


the  above  matrix  equation  takes  the  following  form: 
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(2.33) 
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So  the  above  matrix  equation  gives  the  3  different 
responses  of  the  programm6d)le  filter  as  functions  of  the  Yl, 
.  .  .  Y7  admittances  and  Al,  A2  (Gains)  the  two  Op.  Amps.  In 
ideal  case  [Al=A2=oo]  the  equation  takes  the  form: 
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(2.34) 


which  expressed  the  different  realizations  discussed  by  [31]. 
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Now  we  can  express  the  above  matrix  equation  for  the 
nonideal  case  in  which  A1  and  A2  are  finite  and  frequency 
dependent  where  W1  and  W2  are  the  Gain 

Bandwidth  Products  (GBWP)  of  the  Op.  Amps.  Al,  A2, 
respectively. 


(2.35) 


1.  Low  Pass  (LP)  Realization 


Can  be  obtained  through  T2(s)  by  substituting 
Y6=Y7=0,  and  if  we  consider  the  ideal  case  where  (Al=A2=oo). 


Y3=sC3+G3,  Y4=G4,  Y5=G5,  Y8=G8  then  (2.36)  takes  the  form 


?o  -  I  C^Ci(se 


(2.38) 


<G,  GA  6e 

^3  3  66 


and 


.  6ji 
Qp  (I3 


(2.39) 


(2.40) 


2.  High  Pass  (HP)  Realization 

This  one  is  obtained  through  Tl(s)  substituting  the 
following  values  of  the  admittances;  Y1=G1,  Y2=G2,  Y3s:sC3, 
Yr=G4,  Y5=0,  Y6=G6,  Y7=sC7,  and  Y8=G8,  and  if  we  assume  ideal 
case  (Al=A2->  00),  the 

TlCs:)^  CC■:^^GA')  _  •  (2.41) 


expresses  the  H.P.  filter  of  (2.5)  where 


(si  G<x  Gc. 

G2  C5 


(2.42) 
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Op 


6-3 


(2.44) 


3 .  Band  Pass  (BP)  Realization 

Ttiis  is  also  derived  from  Tl(s)  by  substituting  tlie 
following  values  of  admittances:  yi=Gl,  Y2=G2,  Y3=sC3, 

Y4=G4,  Y5=0,  Y6=G6,  Y7=G7,  and  Y8=sC8,  and  if  we  assume  ideal 
case  (A1=A2 — >oo),  then 


<50  (st,  ^S<SiG  \Cz>  t-Ci  Cq  (j7S^ 

which  is  the  form  of  a  BP  transfer  function  as  given  by  (2.5) 
where 


6-6 
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(2.46) 


Gf 


(2.47) 
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Qj 


(2.48) 


Notch  (N)  Realization 


T2(s)  expressed  the  Notch  response  if  the  following 
substitutions  have  been  made.  Y1=G1,  Y2=G2,  Y3=sC3,  Y4=G4, 
Y6=0,  Y5=G5,  Y7=sC7,  Y8=G8 ,  and  again  assuming  ideal  case 
( A1=A2  —  >oo)  , 


^  ($3 -t- S  v-<S>  <3-+-G5 


(2.49) 


which  is  the  form  of  a  N  transfer  function  as  defined  by 
(2.7)  where 


^3  Cl  Cs2. 


(2.50) 
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(2.51) 
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(2.52) 


5 .  All  Pass  (AP)  Realization 

This  is  derived  from  Tl(s)  with  the  following 
admittances  substitutions.  Y1=G1,  Y2=G2,  Y3=sC3.  Y4=G4. 


'C 


Y5=G5,  Y6=0,  Y7=sC7  and  Y8=G8 ,  and  if  once  more  we  assume 
ideal  case,  then 


^  S^S  Gs  (3-5 

^  ^ ^3(3268  V  G5G1  G4 


(2.53) 


which  is  the  response  of  an  All  Pass  filter  as  (2.8) 


indicates  where: 
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(2.54) 


^i:p  = 
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(2.55) 


Table  2.1  shows  all  the  realizations  proposed  by 
[31],  In  our  research  for  designing  a  programmable  filter 
the  No.  1,  3,  7,  9,  and  12  realizations  were  used  since  they 
offer  the  minimum  admittance  elements  change  to  shift  from 
one  to  another. 

D.  SENSITIVITY  ANALYSIS 

Consider  the  CGIC  circuit  shown  in  Figure  (2.2(a))  [12]. 
Assuming  ideal  Op.  Amps.,  the  chain  matrix  of  the  CGIC  can  be 
obtained  as 


'.•’■y 'o'* 


Fig.  2.3C®0The  CGIC  Implimentation  Using 
Op.  Amps. 

(b)  Symbolic  Representation  of  the 
CGIC  with  created  ports  3G  and  4G 

(c)  The  Basic  Conf irguration. 
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Table  2.1  -  Elements  identification  for 
Realizing  the  Most  Commonly  Used  Transfer  Fuctions 


= 


1  0 
L  ^ 


(2.56) 


where  h(s),  the  admittance  conversion  function,  is  given  by: 
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Two  new  ports  can  be  created  across  terminals  3G  and  4G 
as  shown  symbolically  in  Fig.  (2.3  (b)). 

A  synthesis  procedure  is  now  described  which  uses  the 
configuration  of  Fig.  (2.3(c)).  The  transfer  functions 
between  the  input  and  output  terminals  2,  3,  and  4  are 
readily  obtained  as 


r  Vi-:  1 Dci) 
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(2.58) 


(2.59) 


(2.60) 


‘  V* 


where 


(2.61) 


The  conversion  function  h(s)  and  Y5-Y8  can  be  selected  in 
many  different  ways  and  it  is  found  that  any  second-order 
transfer  function  can  be  realized  [12]. 

Letting 

''!l  ■:  sC  2.  ^  Gx 

where  i=l,  2,  3,  or  4,  we  have  from  (2.57) 

llci) CiiCi  *Gz')CsCi  I  CsCi  '-GO  lsC^  +  GO)  (2.62 ) 

Clearly  by  omitting  one  or  more  conductances  and/or 
capacitances  a  number  of  specific  conversion  functions  can  be 
generated. 

Most  frequently,  filters  are  designed  by  using 
Butterworth,  Chebychev,  Bessel  or  elliptic  approximations  in 
which  the  transmission  zeros  are  located  at  the  origin, 
imaginary  axis  or  at  infinity.  Consequently,  the  transfer 


function  can  be  expressed  as  a  product  of  a  second-order 
transfer  functions  of  the  form 


TCS^  r 


Oi  Oa  s  do 
bj  S'vbi^*-b<p 


(2.63) 


where  al=a2=0,  ao=al=0,  ao=a2=0  or  al=0  for  low  pass  (LP), 
high  pass  (HP),  band  pass  (BP)  or  notch  (N)  section, 
respectively. 

The  coefficients  ai's  of  T(s)  for  these  sections  are  all 

positive.  These  sections  can  be  realized  by  choosing  h(s)  in 

2 

a  simple  manner  such  as  kj^s,  k2S,  k^s  +k^s  or  their 
reciprocals.  The  different  ki's  (i=l,  2,  3,  or  4)  are 
positive  constants.  By  comparing  (2.58-2.61)  with  (2.63), 
circuits  1-10  in  Table  2.1  can  be  obtained.  Circuits  3,  4 
and  7  can  be  regarded  as  realizations  of  simple  RLC  networks 
[31]. 

All  pass  transfer  functions  are  often  needed  for  delay 
equalization  and  these  can  be  realized  by  using  second-order 
transfer  functions  of  the  type  give  by  (2.63)  where  a2=b2, 
al=bl,  and  ao-bo.  Second-order  sections  of  this  class  can  be 
obtained  from  circuits  11  and  12  of  Table  2.1. 

Figures  (2.3(a)),  (2.3(b)),  and  Table  2.1  show  that  with 
the  exception  of  circuit  10,  the  response  is  obtained  from 
the  output  of  an  operational  amplifier.  Owing  to  the  low 
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output  resistance  of  the  amplifier,  any  number  of  sections 
can  be  cascaded  without  isolating  amplifiers. 

An  important  criterion  of  a  realization  is  its 
sensitivity  to  element  variations.  The  pole  Q  factor  and  the 
undamped  frequency  of  oscillation  from  the  transfer  function 
of  (2.63)  are  defined  as 


r  (bT^  /bi 


)  {'^1^2 


(2.64) 


For  a  Notch  section,  the  Notch  frequency  is  defined  by 


i  (.Qcl  Q'j') 


and  the  multiplier  constant  can  be  taken  to  be 


HmzOo  be  or  Oi  I  bi. 


for  tjJn>oip  or  uc>n<obp,  respectively. 

Similarly  for  the  LP,  HP,  and  BP  sections 
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(2.65) 


(2.66) 


(2.67) 


For  an  All  Pass  Section,  let 


^  lC3=(^0©|Qa) '^o„4  \\A,?-03|b2 


(2.68) 


The  sensitivity  of  a  quantity  x  with  respect  to 
variations  in  an  element  e  is  given  by 


For  ideal  amplifiers,  the  use  of  (2. 64) -(2. 68)  and  Table 
(2.1)  leads  to 


(2.69) 


where  x  represents  any  one  of  the  quantities  defined  by 
( 2.64)-( 2.68)  and  e  represents  any  capacitance  of 
conductance.  In  addition,  it  can  be  shown  that 


For  amplifiers  with  a  finite  open- loop  gain  A,  according 
to  [31]  the  circuit  of  (Fig.  1(c))  gives 


N  v.s/' 


(2.71) 


where  k=2,  3,  4,  and 

rc>:)=Fi.  yv  tfs '(3  I- Cl  *■  I  (2.72) 
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Consider  realizations  in  which  h(s)  =  kls,  such  as  the 
circuits  7(BP),  3(HP),  9,  10(N),  where 

(5:2 1  tb  I S  Xd- (s*. Nt-(sfe  I  h- SC?  vCs’  Odil 
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For  real  amplifier  gains  such  that 

Qvii  Ao /)  L 


(.T  17)  gives 
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From  (2.64),  (2.73),  and  (2.74),  trie  Q-f actor  and  the 

undaunped  frequency  of  oscillation  can  be  obtained  as 


'Ij  ’ll  . 

^'Ao  J  i  /  I  ~C3  >>2  Ao  \ 


1 


i  (Ci4Ce)Ao'. 


(2.75) 


where 
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(2.76) 

'^5-  (C^  , 

^5  =  ood  'si^ 

The  sensitivities  of  Qpa  and  pa  with  respect  to  the 
amplification  Ao  can  be  written  as 
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The  use  of  (2.77)  and  (2.78)  leads  to 


(2.78) 


by  assuming  that 


Eq.(2.78)  reduces  to 
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(2.79) 


(2.80) 


Straight  forward  differentiation  shows  that  S*^  is  minimum 
when 
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From  (2.79)  and  (2.81)  the  analysis  is  valid  provided 


that  4  G'^  *4- 
practice. 


which  is  clearly  satisfied  in 


a 


From  (2.80)  and  (2.81)  the  minimum  sensitivity  to 
variations  in  Ao  is  derived  as 


/Ac 


(2.82) 
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The  corresponding  value  of  is  given  by 

Ao 

..  1 


Ao  =  ~ 


Ao^p 


LP  realizations  as  circuits  1  and  2  can  be  obtained  by  using 
a  conversion  function  of  the  form  or  its  reciprocals.  The 
admittances  Y1  to  Y8  are  chosen  as 

?Z2,  O  .Vf:.  (2.84) 


in  order  to  obtain  a  conversion  function  of  the  form  \ 


NS  ■- 


Now  Qpa  and  (vpa  are  obtained  as 


1  ^  ^  ;  . . .  . - 


(2.85) 
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(2.86) 
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The  sensitivity  of  Qpa  with  respect  to  variations  in  Ao 
can  be  minimized  following  the  approach  used  earlier.  It  is 


found  that  for  minimum  sensitivitv[ 311 


(ja  -  (^B 


^  65 


(2.88) 


The  minimum  value  of  S, 


Qpa 


Ao 


can  be  shown  to  be 


(2.89) 


and  the  corresponding  value  of  is  given  by 


c;'^?cc 
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(2.90) 


The  above  sensitivity  analysis  can  be  extended  to 
realizations  using  any  other  type  of  conversion  function 
C31]. 

Equation  (2.69)  shows  that  the  sensitivities  to  passive 
element  variations  are  independent  of  the  selectivity. 
Furthermore,  the  sensitivities  with  respect  to  variations  in 
the  amplifier  gain  are  low.  The  proposed  realizations  are 


seen  to  have  similar  sensitivity  properties  as  the  low 
sensitivity  realizations  reported  in  [21-27]. 

E.  STABILITY 


It  has  been  shown  elsewhere  [20]  that  some  networks  using 
GIC's  can  be  conditionally  stable  where  a  circuit  can  lock  in 
an  unstable  mode  during  activation  (just  after  switching  on 
the  power  supply).  In  this  section  the  stability  properties 
of  the  configuration  show  in  Fig.  (2.2(c))  are  exeunined. 

The  natural  frequencies  of  the  circuit  in  Fig. (2.3(c)) 
are  the  zeros  of  the  characteristic  polynomial  D9s)  as  given 
by  (2.72).  The  differential  open-loop  gain  of  a  frequency 
compensated  Op.  Amp.,  in  a  bounded  frequency  range  C=c. £  ;  v  : 

A-i, 

where  Ao  and  wc  are  the  d.c.  gain  and  the  cutoff  frequency 
respectively,  and  i  ;\o  4  Ai'vo;,  .  In  the  frequency  range 
i.<<-£c  the  cimplifier  gains  Al  and  A2  can  be  assvuned  to  be 
real.  For  any  second-order  transfer  function  the  coefficients 
of  D(s)  are  seen  to  remain  positive  for  any  attainable  pair 
of  Al  and  A2.  this  is  due  to  the  absence  of  negative  terms 
in  D(s).  Therefore,  the  zeros  of  D(s)  will  remain  the 
left-half  s-plane  and  low  frequency  unstable  modes  cannot 
arise  during  activation. 
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III.  PROGRAMMABLE  GIC  FILTER 

A.  GENERAL 

Signal  processing  devices  evolved  considerably  over  the 
last  several  years.  The  progress  was  motivated  by  the 
advancement  in  film  and  semiconductor  technologies,  as  well 
as  the  continuous  upgrading  of  systems  specifications  to  take 
advantage  of  the  available  technologies  to  the  limits. 

Linear  filtering  finds  many  applications,  such  as  speech 
processing  (recognition  or  synthesis),  geology, 
instrumentation,  communications,  process  control,  adapting 
balancing,  etc.  There  has  been  much  emphasis  on  performing 
the  filter  function  digitally,  largely  because  of  the  ease  of 
varying  and  optimizing  the  transfer.  However,  and  for  many 
reasons,  such  as  cast  size  signal  processing  complexity,  and 
bandwidth,  it  would  be  desirable  to  perform  the  filter 
function  with  linear  components,  yet  retain  the  flexibility 
of  varying  the  filter  parameters  digitally. 

Recently,  several  advantages  of  combining  linear 
components  (amplifiers  and  capacitors)  and  nonlinear  elements 
(switches)  have  been  demonstrated  using  MOS  switched 
capacitor  techniques  [31,  32].  Here,  we  are  presenting  the 
results  of  realizing  a  continuous  active  device  using  linear 
elements  and  switches  controlled  by  digital  signals  to 
achieve  fully  programmable  filters. 


Our  research  addressed  two  different  aspects  of 
programmability  namely; 

(1)  Programming  the  filter  topology  using  a  minimal  set  of 
elements  to  obtain  any  type  of  filtering  function  desired, 
e.g.,  LP,  HP,  BP,  N  and  AP. 

(2)  Programming  the  filter's  transfer  function  parameters, 
(pole  resonant  frequency  iop  and  quality  factor  Qp)  for  a 
chosen  type  of  filtering  function. 

B.  THE  PROPOSED  GIG  PROGRAMMABLE  FILTER 

The  basic  active  network  considered  as  the  heart  of  the 
GIG  programmable  filter  is  the  GIG  structure  [31]  of 
Fig. (3.1),  whose  superior  performance  was  established  in  the 
literature  [10,32].  The  filter  transfer  function  was  derived 
using  loop  analysis  in  Ghapter  II, 

Table  (3.1),  illustrates  that  for  any  of  the  LP,  HP,  BP, 
N  and  AP  realizations,  five  resistors,  two  capacitors  and  two 
Op.  Amps.,  are  required,  also,  the  transfer  function  of  each 
realization  is  shown.  The  passive  elements  are  connected  to 
the  different  nodes,  shown  in  fig. (3.2),  for  the  different 
realizations.  A  set  of  MOS  bilateral  switches  controlled  by 
a  digital  binary  word,  are  used  to  interchange  the  elements 
to  achieve  the  different  types  of  filter  realization  shown  in 
Fig. (3.3).  The  truth  table  of  the  switch  control  logics  is 
shown  in  Table  (3.2).  Fig. (3.3(a))  illustrates  the  GMOS 
logical  circuit  for  realizing  this  truth  table.  While  four 
of  the  resistors  are  equal  and  of  value  R  each,  the  fifth 
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The  CMOS  Logic  Diagram  Used  to  Control  the 
Analog  CMOS  Switches  of  Fig.  and  to 
Realize  the  Truth  Table 


resistor  is  the  Qp  determining  resistor  and  of  value 
Rq=RQp.  The  two  capacitors  are  equal  and  of  value  C=l/WpR. 
each.  The  two  equal  banks  of  capacitors  are  used  to  control 


up.  Each  bank  contains  n  binary  weighted  capacitors 
connected  in  series  through  analog  CMOS  switches  as  shown  in 
Fig. (3.5).  Using  a  digital  binary  word  of  n  bits  to  control 
Wp,  2^  different  values  of  C  will  result  at  the  2  terminals 
of  both  capacitors  banks  that  correspond  to  2^  different 
values  of  p.  Using  a  similar  technique  the  value  of  Rq  can 
be  controlled  through  a  bank  of  m  binary  weighted  resistors 
in  series,  though  analog  CMOS  switches  as  shown  in  Fig. (3.6). 
Using  a  digital  binary  word  of  m  bits  to  control  Qp,  2*^ 
different  values  of  Rq  can  be  achieved  that  correspond  to  2*^ 
different  values  of  Qp.  Thus,  full  independent  control  of 
the  pole  pair  (Op  and  Qp  are  achieved  by  programming  the 
switches  to  obtain  the  corresponding  C  and  Rp.  It  can  be 
easily  shown  that  with  minor  modifications,  an  additional 
programmable  element  can  be  added  for  the  control  of  the 
notch  frequency. 

C.  THE  REALIZED  GIC  PROGRAMMABLE  FILTER 


A  complete  circuit  diagram  of  the  constructed  GIC  filter 
is  shown  in  Fig.  (3.4).  The  values  of  m  and  n  were  selected 
to  m=n=4.  Thus,  15  different  values  of  lOp  (fp)  and  Qp  were 
obtained  as  it  is  illustrated  at  the  corresponding  Table 
(3.3)  and  (3.4).  the  ( designed )  banks  of  the  resistors  for 
the  control  of  Qp  and  the  capacitors  for  the  control  of  oDp 
(fp)  along  with  their  control  switches  are  shown 
correspondingly  in  Fig. (3.7)  and  Fig. (3.8). 
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Fig.  3.8  -  The  Constructed  Block  of 
Resistance  Controlling  O 


IV.  COMPUTER  SIMULATION  OF  THE  PROGRAMMABLE  FILTER 


A.  INTRODUCTION 

In  order  to  observe  the  theoretical  responses  of  the 
different  realizations  as  well  as  to  compare  them  with  the 
experimental  measurements,  two  computer  programs  were  written 
in  Fortran,  those  programs  are  shown  in  Appendices  A  and  B. 
The  first  one  simulates  the  progremmiable  filter's  responses 
Tl(s),  T2(s)  and  T3(s)  as  functions  of  the  network's 
admittances  Yi,  i=l  .  .  .8,  and  the  "constant"  Op. Amps,  gains 
A1  and  A2.  A  realistic  value  of  A2=A2=10^  was  given  to  the 
above  gains  corresponding  to  741  Op.  Amps,  used  in  the 
experiment.  the  second  program  simulates  the  programmable 
filter's  responses  Tl(s),  T2(s)  and  T3(s)  as  functions  of  the 
network's  admittances  Yi,  i=l  ...  8  and  the  frequency 

dependent  Op.  Amps,  gains  A1  and  A2.  In  this  second  case  a 
single  pole  approximation  value  of  wils  was  assigned  for 
Ai(i=l,  2).  Where,  Wi  is  the  GBWP  of  the  Op.  Amps,  used,  a 

g 

value  of  W1=W2=2'X  xlO  was  given  to  the  above  Wi(i=l,  2) 
corresponding  to  741  Op.  Amps,  used  throughout  the  research. 

To  be  able  to  compare  the  computer  simulation  results 
with  the  experimental  ones  obtained,  the  same  values  of 
admittances  were  used  as  input  data.  That  means  that  the 
same  values  of  R,  Rq  and  C  according  to  Chapters  III  and  V 
were  used.  With  the  value  of  R  selected  to  be  1.6k  Table 
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(4.1)  illustrates  the  different  values  of  Rgs  used  in  the 
simulation  to  realize  the  different  values  of  Qp  and  Table 

(4.2)  illustrates  the  different  values  of  capacitors  used  in 
the  simulation  to  realize  the  different  frequencies  Wps. 

The  "DO  CASE  I"  command  of  Fortran  simulated  the  digital 
logic  (including  the  control  switches)  used  to  realize  the 
different  responses  of  the  programmable  filter  that  is  LP, 
HP,  BP,  N  and  AP.  The  frequency's  (i;ap)  translations  and  the 
different  QP  values  were  simulated  by  changing  the  values  of 
Rq  and  Cs  at  every  run  of  the  program. 

The  above  programs  also  simulated  the  transfer  function 
of  two  cascaded  QIC  programmable  filters  as  it  will  be 
discussed  later  in  Chapter  VI.  Each  of  the  filters  could 
have  been  at  different  realization  (as  well  as  at  different 
p  and  Qp)  relative  to  the  other  one.  The  25  possible 
combinations  of  transfer  function  realizations  are  shown  in 
Table  4.3. 

B.  SIMULATION  RESPONSE(S) 

1 .  Low  Pass  (LP)  Realization 

Using  the  elements  values  prescribed  in  Table  (3.1) 
yields  the  following  L.P.  transfer  function; 
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The  25  Possible  Combinations  of 
Transfer  Function  Realizations 


Fig.  4.1  illustrates  the  theoretical  "ideal"  LPF 
magnitude  response  for  fp=3.8Khz  and  for  3  different  values 
of  Qp.  The  simulation  results  match  the  equations  of  (4.2). 

Fig.  4.2  illustrates  the  theoretical  "ideal"  LPF 
magnitude  response  for  Qp=2  and  for  3  different  values  ofu,p. 
the  Op.  Amps,  gains  Ai,  (i=l,  2)  are  frequency  depended. 
This  dependence  affects  the  magnitude  of  the  filter  and 
causes  a  frequency  shift  from  the  theoretical  value  of  the 
ideal's  case.  Figs. (4.3)  and  (4.4)  illustrates  the  ideal  vs. 
nonideal  theoretical  LPF  amplitude  responses. 

Data  extracted  from  Figs. (4.3)  and  (4.4)  are 
illustrated  in  Tables  (4.4(a))  and  (4.4(b))  simultaneously. 


8a  a 


LPJ  AMPLRESPONSE  (Q=2) 


-  ideal 

_  nonideal 

Fig.  4.3  Ideal  .vs.  nonideal  L.F.F.  amplitude  response  for  Q=2 
and  frequencies  (1.99K,  7.96K,  30.16K  and  117.lK)h> 
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Data  From  Fig.  (4.4)  LPF  Characteristic 
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the  frequency  increases)  while  Figs.  (4.8)  and  (4.9) 
illustrate  the  frequency  shift  of  the  amplitude  response  due 
to  that  dependency. 

From  the  date  of  Fig.  (4.8),  Table  (4.5)  were 

constructed. 

3 .  Band  Pass  (BP)  Realization 

Using  element  values  from  Table  3.1  the  following  BP 
transfer  function  is  achieved: 


takes  the  following  value  at  top. 


(4.7) 


T-llitOO  -  0. 
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which  has  constant  magnitude  of  6dB 


(4.9) 
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Fig.  (4.10)  illustrates  the  theoretical  BPF  magnitude 
response  for  fp=3.83  KHZ  for  different  values  of  Qp.  These 
agree  with  the  (4.9)  equation  since  as  it  is  indicated  by  the 
simulation  plot  the  amplitude  is  constant  and  independent  of 
Qp. 

Fig.  (4.11)  illustrates  the  above  concept  but  at 
fp=15.1KHZ. 

Fig.  (4.12)  and  (4.13)  illustrate  the  theoretical 
"ideal”  BPF  amplitude  response  for  Qp=2  and  for  different 
values  of  p.  As  it  is  indicated  from  Fig.  (4.13)  the 
amplitude  remain  constant  even  at  very  high  frequencies  (10 
HZ).  But  with  Ai,  (i=l,  2)  depending  on  frequency  the 
amplitude  decreases  as  the  frequency  increases.  This  is 
indicated  in  Figs.  (4.14)  and  (4.15)  which  describe  the  BPF 
amplitude  response  plots  for  Qa4  and  Q=1  respectively  and  for 
different  frequencies.  The  frequency  dependence  of  Ai,  (i=l, 
2)  creates  a  frequency  shift  from  the  ideal  theoretical  value 
which  is  indicated  in  Figs.  (4.16)  and  (4.17).  Table  (4.6) 
illustrates  the  data  extracted  from  figs.  (4.16)  and  (4.17). 

4.  Notch  (N)  Realization 

Using  the  admittances  value  of  Table  (3.1),  the 


following  Notch  transfer  function  can  be  achieved: 


B.P.F  AMPLITUDE  RESP0NSE(Q=4) 


Fig.  4.12  -  'Ideal'  BPF  Amplitude  Response. 
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Fig.  4.13 


'Ideal'  BPF  Amplitude  Response. 
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where  n  is  the  Notch  frequency. 

At  P  the  transfer  function  takes  the  value 

TiCjooy)  -  (4.11) 

Qp 

Fig. (4.18)  illustrates  the  "ideal"  theoretical  Notch 
filter  amplitude  response  for  a  variety  of  frequencies  and 
constant  Qp(Q=2) ,  while  Fig. (4.19)  for  a  variety  of  Qps  and 
for  constant  p  (fp=3.83KHZ) . 

Fig.  (4.20)  illustrates  the  effect  of  the  frequency 
dependency  of  Ai,  (i=l,  2)  on  the  response. 

5.  All  Pass  (AP)  Realization 

As  proposed  in  Table  (3.1)  using  the  same  admittances 
values  as  Notch.  An  All  Pass  transfer  function  can  be 


derived  as 


which  takes  the  following  values 
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Fig.  4.22 


Ideal  vs.  Nonideal  Phase  Response. 
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V.  REALIZATION  OF  PROGRAMMABLE  GIC  FILTER 


A.  EXPERIMENTAL  RESULTS 

After  the  circuit  of  Figs.  (3.10),  (3.9),  (3.8),  (3.7), 

was  constructed  a  variety  of  measurements  were  taken  in  order 
to  study  the  response  of  the  network  to  the  different  inputs 
(control  bitwords).  To  observe  the  affect  of  the  control 
switches  which  introduce  a  resistance  of  80  each  at  CLOSED 
position,  two  values  of  Rs,  Rq,  and  Cs  were  used  with  one 
decade  difference  in  magnitude.  That  means  that  R  was  given 
the  values  of  1.6K  (as  discussed  in  Chapter  III)  and  16R,  the 
four  resistors  that  consisted  the  Rq  bank  were  of  values 
(1.6K,  3.2K,  6.4K,  12. 8K)  in  the  first  case  and  (16K,  32K, 
64K  and  128K)  in  the  second  one,  and  that  the  capacitor 
bank's  capacitors  were  chosen  of  values  (lOOnF,  50nF,  12, 

5nF,  llnF)  and  (lOnF,  5nF,  1.2nF,  O.lnF)  accordingly,  to  be 
able  to  keep  the  range  of  frequencies  as  much  the  same  as 
possible  for  both  cases. 

1.  Low  Pass  Filter 

With  the  topology-control  bit  word  000  the  network 
realized  a  LPF  response.  Fig.  5.1(a)  illustrates  the  LP 
response  for  R=1 . 6Ki  to  a  variety  of  frequencies  for  0=5 
while  5.1(b)  illustrates  the  same  but  for  Q=2.  It  can  be 
observed  from  5.1(a)  (Q=5)  that  the  amplitude  response 


=  0.99khz  (b) 

=  1.98khz 
=  3.83khz 
=  7.96khz 
=  12.0khz 


Fig.  5.1  -  (a)  LPF  Response  for  R=1.6k  (Q=5,  Rq=8.0k) 
(b)  LPF  Response  for  R=1.6k  (Q=2,  Rq=3.2k) 


decreases  while  the  frequency  increases  as  it  was  expected 


from  computer  simulation.  This  is  occurred  up  to  the 
frequency  of  6KHZ;  then  it  started  increasing  with  the 
frequency,  while  at  Fig.  (5.1(b))  (Q=2)  it  remained  constant. 
Fig.  (5.2)  illustrates  (for  R=16KHZ)  that  the  observation  in 
Fig  (5.1)  is  not  any  more  the  case  (for  that  frequency  range) 
and  the  networ]c  responses  the  scune  as  in  computer  simulation 
while  in  Fig.  (5.3)  (Q=5)  the  above  can  be  noticed  again, 
this  is  due  to  the  interference  of  the  control  switches  as  it 
was  discussed  previously.  Figs.  (5.4)  and  (5.5)  illustrate  a 
variation  of  Q  values  for  f=9KHZ  for  the  two  cases  (R=1.6k 
and  R=16K  )  respectively.  A  difference  in  magnitude  can  be 
observed  due  to  the  interference  of  the  control  switches. 
Figs.  (5.6)  and  (5.7)  illustrate  the  same  but  for  F=12.8KHZ. 
Figs.  (5.8)  and  (5.9)  illustrate  the  phase  and  amplitude 
response  of  the  LPf. 

2.  High  Pass  (HP)  Realization 


With  the  topology-control  bitword  001,  the  network 
realized  a  high  pass  filter.  Figs.  (5.10)  and  (5.11) 
illustrate  the  HPF  amplitude  response  for  a  variety  of 
frequencies  and  (Q=2).  It  can  be  observed  (as  in  LPF 
realization)  that  for  this  frequency  range  and  for  R=1.6I& 
the  amplitude  starts  increasing  as  the  frequency  increases, 
decreases  then  again  while  for  R=16KX  this  does  not  occur, 
the  same  observation  submerges  comparing  Fig. (5.14)  and 
(5.15).  Figs.  (5.12)  and  (5.13)  are  the  plots  obtained  for 


2ibl<ik  : 


0 . 99khz 
1.99khz 
3 . 83khz 
7 . 90khz 
12 • 8khz 


Fig.  5.2  - 


LPF  Amplitude  Response 
R=1.6k  (Q=2,  Rq=32k) 


=  0.99khz 
=  1.99khz 
=  9,20khz 
=  7.96khz 
=  IS.Okhz 


Fig.  5.3  -  LPF  Amplitude  Response 
R=16k  (Q=5,  Rq=50k) 
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1  -  f  =  0.99khz 

2  -  f  =  1.99khz 
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4  -  f  =  7.96khz 


f=6.65KHZ  and  a  variety  of  Qs  for  R=1.6K£  and  R=16K2  , 
respectively.  a  difference  of  about  12dB  in  amplitude 
appears.  Finally,  figs  (5.16),  (5.17),  (5.18),  (5.19) 
illustrate  the  phase  and  amplitude  responses  of  the  HP 
realization. 

3.  Band  Pass  (BP)  Realization 

The  010  topology-control  bitword  realizes  the  Band 
Pass  Filter.  Figs.  (5.20)  and  (5.21)  illustrates  the 
amplitude  response  for  set  of  frequencies  and  Q=10.  At  both 
cases  the  amplitude  decreases  while  the  frequency  increases 
until  the  value  of  9KHZ.  Then  starts  increasing  again.  It 
can  be  also  observed  that  for  R=1.6K  the  frequency  deviation 
from  the  theoretical  fp  is  larger. 

Figs.  (5.22)  and  (5.23)  illustrate  the  amplitude 
response  again,  but  for  Q=l.  This  time  the  deviations  from 
the  theoretical  response  (both  of  amplitude  fluctuations  and 
frequency  shift)  are  less. 

Figs.  (5.24),  (5.25),  (5.26),  and  (5.27)  are  also 
plots  of  the  amplitude  response  but  for  variation  of  Q.  A 
difference  of  approximately  IdB  appears  for  the  lowest 
frequency  (3.8KHZ)  and  of  0.5dB  for  the  higher  (9KHZ), 
Finally,  Figs.  ( 5. 28) -(5. 32)  illustrate  the  phase  and 
amplitude  of  the  BPF  response. 

4.  Notch  (N)  Filter  Realization 


With  the  topology-control  bitword  Oil,  a  Notch  filter 
realization  can  be  achieved.  Figs.  (5.33)  and  (5.34) 
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Fig.  5.19  -  HPF  Amplitude  Response  (0=2) 
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Fig.  5.20  -  BPF  Amplitude  Response 
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Fig.  5.21  -  BPF  Amplitude  Response 
(R=1.6k),  Q=10 
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Fig.  5.32  -  BPF  Phase/Amplitude  Response 
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Fig.  5.34  -  Notch  Amplitude  Response 
(R=1.6k) 


illustrate  the  amplitude  response  for  different  frequencies 
and  for  Q=4,  while  Figs.  (5.35)  and  (5.36)  illustrate  the 
amplitude  response  for  f=lKHZ  and  a  variety  of  Qs.  Fig. 
(5.37)  illustrate  the  phase  response  in  addition  to  the 
amplitude  one. 

5 .  All  Pass  (AP)  Realization 

The  topology-control  bitword  100  realizes  the 
All-Pass  filter.  Figs  (5.38),  (5.39)  and  (5.40)  illustrate 
the  plotted  amplitude  and  phase  responses. 

B.  CONCLUSION 

The  constructed  circuit  performed  as  predicted  by  the 
theoretical  analysis  and  the  computer  simulations.  This 
means  that  it  realized  all  the  desired  filtering  transfer 
functions  LP,  HP,  BP,  N,  and  AP.  The  effect  of  interference 
of  the  control  switches  nonideal  performances  which  is  more 
severe  at  high  Qps  can  be  minimized  by  increasing  the  values 
of  R  (meaning  that  the  values  of  Rg  also  increases  and  the 


values  of  Cs  decreases ) . 


2.  Amplitude  Response. 


KbPKUUUUk. 


ISS&il^S&SSSSSSSESSSl 


m 


HtPKOUU»_ 


MFR  FCTH  MKP^  OFF  SCFILE 
XFR  FC T  H  =  0  o  CENTER 


dB  0  I 
0  o  .  D  I  'J 


By  cascading  two  or  more  programmable  filters,  higher 
order  transfer  functions  can  be  obtained.  Fig.  (6.1) 


illustrates  two  cascaded  GIC  programmable  filters, 


Fig.  (6.1)  Two  Cascaded  GIC  Programmable  Filters. 

(Each  Black  box  Stands  for  the  Network 
of  Fig.  (3.10) 


25  different  combinations  of  the  two  individual  transfer 
functions  Ti(s)  and  Tj(s)  can  obtained  as  it  is  indicated  in 
Table  (4.3). 

For  Lp-LP  combination  and  for  ideal  theoretical  case 
(A1=A2 — >oo)  a  fourth  order  low  pass  filter  can  be  obtained 
with  transfer  function; 


The  computer  simulation  of  the  fourth  order  transfer 
function  of  a  nonideal  theoretical  filter  vs.  the  second 


order  case  is  illustrated  at  Fig.  (6.2).  The  experimental 
results  taken  from  the  properly  designed  and  built  circuit  as 
illustrated  in  Fig.  (6.3)  are  shown  in  Fig.  (6.4). 

Using  the  same  procedure  as  above  a  fourth  order  HP-HP 
combination  for  the  ideal  theoretical  response  is  given  by 


\\Ci)  ^  4  ^ _ 

which  at  Wp  takes  the  complex  value  of 


(6.2) 


(6.3) 


with  magnitude  of 

|Ti(j0.ip)j=  40^^  =  40tog(70^}<)l,  (6.4) 


Fig.  (6.5)  illustrates  the  fourth  order  nonideal 
theoretical  HP  filter  response  vs.  the  second  order  one. 
Both  at  Qp=2  and  fp=8KHZ ,  while  fig.  (6.6)  illustrates  the 
experimental  responses.  For  Qp=2  according  to  (6.4)  and 


LP.F  AMPL  RESP0NSE(Q=^/F=4K) 


Fig.  6.3  -  Network  with  Logic  for  Realizing 

np  to  8th  Order  LP,  HP,  BP,  N,  and 
AP  Transfer  Functions. 


(2.4)  a  difference  of  12dB  were  expected  between  the  second 
and  fourth  order  filter  which  is  approximately  the  case  in 
both  computer  simulation  and  experimental  results. 

The  BP-Bp  combination  leads  to  the  theoretical  "ideal" 
fourth  order  transfer  function: 


(6.5) 


which  takes  the  following  value  at  a'p: 

h  (j  (6.6) 

According  to  this,  a  magnitude  response  of  12  dB 
approximately  at  (Cp,  was  expected  from  both  experimental  and 
computer  simulation  results.  Fig.  (6.7)  illustrates  the 
simulated  response  of  a  fourth  order  nonideal  HP  filter  vs.  a 
second-order  one.  A  difference  of  5.8  dB  instead  of  6  dB  can 
be  observed.  Fig.  (6.8)  Illustrates  the  experimental 
response,  where  a  difference  of  5.4  dB  appears  basically  due 
to  the  interference  of  the  control  switches  and*  the 
nonideally  matched  values  of  the  capacitors  which  control  the 
fp  selection. 
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B.P.F  AMPU  RESP0NSE(Q=4/P=8K) 


The  N-N  combination  for  the  ideal  theoretical  case 
(A1=A2 — >oo)  results  the  following  transfer  function: 


T2C5.')=  ( 


(6.7) 


where  (On  is  the  Notch  frequency. 

Fig.  (6.9)  illustrates  the  fourth  order  "nonideal"  Notch 
filter  simulated  response  vs.  the  second  order  one. 

The  AP“Ap  combination  for  the  theoretical  ideal  case 
result  in  the  following  transfer  function: 


TaCs)  = 


(6.8) 


which  talces  the  following  values  for  S=0  and  oo: 


TiLio'lr  1 


with  amplitude  and  phase  of 


(6.9) 


TiLiCD^  -  a.i5 


(6.11) 


with  amplitude  and  phase  of 

1T(  i  =  7<lb  ,  LJlk2^  =  0'“ib  (6.12) 


The  computer  simulation  of  nonideal  fourth  order  all  pass 
filter  vs.  a  second  order  one  which  is  illustrating  at  Fig. 
(6.10)  matches  the  above. 

Fig.  (6.11)  illustrates  the  resulting  Chebychev  filter 
from  a  BP-BP  combination  with  Qp=4  and  different  frequencies 
(fpl=8KHZ  fp2=10KHZ),  while  Figs.  (6.12)  and  (6.13) 
illustrate  the  resulting  Chebychev  filters  from  the  designed 
and  built  circuit.  Fig.  (6.14)  illustrates  the  resulting 
response  from  a  HP-LP  combination. 
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Fig.  6.10  -  APF  Fourth  Order  vs.  Second  Order 
Amplitude  Response  from  Computer 
Simulation. 
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A.  BACKGROUND 


1.  General  Description  of  Frequency  Hopping  Signals 

Frequency  hopping  is  a  spread  spectrum  modulation 
technique  used  to  generate  many  possible  carrier  frequencies 
over  a  large  bandwidth.  Of  all  the  possible  carrier 
frequencies,  only  one  is  selected  at  a  given  time.  However, 
all  frequencies  are  eventually  selected  during  some  time 
interval . 


Frequency  hopping  (FH)  may  be  pictured  as  an  RF 
carrier  whose  center  frequency  is  "hopped”  over  many 
frequencies.  The  hopping  may  be  either  in  a  simple  sequence 
or  a  pseudorandom  sequence. 

The  hopping  rate  of  a  frequency  hopping  system  does 
not  affect  the  bandwidth  of  the  output  spectrum.  In  a  direct 
sequence  system  the  chip  rate  determines  the  total  bandwidth. 
In  a  frequency  hopping  system,  however,  the  bandwidth  is 
determined  by  the  highest  and  lowest  frequencies  of  the 
frequency  hopped  carriers.  For  example,  if  the  highest 
frequency  carrier  is  at  15  MHz  and  the  lowest  frequency 
carrier  is  at  10  MHz,  the  total  signal  bandwidth  is  5  MHz. 
This  is  the  bandwidth  regardless  of  the  hopping  rate.  This 
allows  wideband  spread  spectrvun  signal  generation  at  low 
hopping  rates. 


Signal  Generation 


Frequency-hopped  signals  may  be  generated  in  several 
ways.  The  different  methods  are  classified  into  two  groups: 

(1)  Direct  synthesis,  and 

(2)  Indirect  synthesis. 

One  important  aspect  of  frequency  hopping  synthesis 
is  coherency.  coherent  signal  synthesis  is  defined  as  the 
establishment  of  a  known  and  repeatable  phase  each  time  a  new 
frequency  is  hopped  to.  Non-coherent  signal  synthesis  is 
defined  as  the  establishment  of  a  random  or  unknown  phase 
each  time  a  new  frequency  is  hopped  to.  some  techniques, 
direct  or  indirect,  of  signal  generation  can  be  used  as  a 
coherent  frequency  source.  In  other  techniques,  the  changing 
of  frequencies  creates  non-coherent  sources. 

If  a  frequency  hop  system  is  a  coherent,  it  will  have 
a  signal- to-noise  advantage  over  a  non-coherent  system. 

a.  Direct  Synthesis 

The  direct  approach  to  signal  synthesis  utilizes 
techniques  which  enable  direct  synthesis  of  different 
frequencies.  Examples  of  direct  synthesis  techniques  are: 

(1)  Frequency  mixing,  and 

(2)  Surface  acoustic  wave  devices. 

Frequency  mixing  for  single  synthesis  is  a  common 
technique  used  to  generate  many  different  frequencies.  An 
example  of  the  frequency  mixing  technique  is  show  in  Fig. 


Fimr.  7 . 1  BLOCK  DI/\CR/W1  OF  A  FREQUENCY  HOP  MODEM 

A  pteudortndofli  cod«  9«n«rPtor  tvlectl  one  of  niny  trantmU  frcqucncici  durfnq  a  imall 
tliM  Interval.  Il<e  traffic  modulates  the  carrier  frequency  which  Is  spread  over  many 
different  frequencies  hy  the  hopping  action,  the  receiver  dehops  the  input  signal 
Into  a  narrowband  If.  The  code  synchronizer  locks  the  pseudorandom  code  generator  In 
the  receiver  to  the  received  signal.  A  data  demodulator  removes  the  traffic  from  the 
If  ami'lifier  output. 
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In  Fig.  (7.2),  an  RF  switch  selects  one  of 
several  frequency  inputs.  Two  of  these  input  signals  of 
different  frequencies  are  multiplied  together  to  generate  a 
new  output  frequency.  The  device  used  to  multiply  the  two 
signals  together  is  called  a  frequency  mixer.  when  two 
frequencies  are  mixed,  the  sum  and  the  difference  of  the 
frequencies  are  generated.  In  order  to  select  only  one  of 
these  frequencies,  a  "filter  is  used  to  reject  the  unwanted 
frequency."  A  filter  tuned  to  the  desired  frequency  would 
allow  selection  of  that  frequency  while  rejecting  the  other. 
By  selecting  the  mixing  frequencies  in  the  proper  order,  the 
output  frequency  can  be  stepped  through  several  different 
frequencies.  At  each  frequency  mixer  output,  a  filters  is 
required  to  reject  unwanted  frequencies.  The  filters  may 
require  a  short  time  for  the  signal  to  stabilize  after  it  is 
selected.  The  time  required  for  the  filter  to  stabilize  at 
each  new  frequency  may  ultimately  determine  the  maximum 
hopping  rate  of  the  direct  frequency  synthesizer, 
b.  Indirect  Synthesis 

The  indirect  method  of  signal  synthesis  is 
defined  as  frequency  synthesis  through  the  use  of 
phase-locked  oscillators.  One  common  indirect  method  for 
synthesis  is  shown  in  Fig.  (7.3). 


In  this  circuit  a  phase-locked  loop  is  used  to 
generate  the  numerous  carrier  frequencies.  The  phase-locked 
loop  has  an  internal  oscillator  whose  output  frequency,  Fo, 
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FIGURE  7.3  INDIRECT  SIGNAL  SYNTHESIS  USING  A  PHASE-LOCKED  LOOP 

When  the  reference  frequency,  Is  at  the  same  frequency  as  F,*,  the  voltage- 
controlled  oscillator  produces!  constant  output  frequency,  Fq.'  Since  F»  Is 
divided  by  .the  programmable  divider,  Tq  Is  equal  to  N  times  F, .  8y*chang1ng 
the  divider  ratio,  N,  many  output  frequencies  are  possible. 


is  shovm  in  Figure  (7.3)  as  the  phase- locked  loop  output.  The 
divide-by-N  circuit  divides  this  oscillator  frequency  by  a 
selected  number,  N. 

The  phase- locked  loop  internally  adjust  its 
output  frequency  Fo  so  that  FI*  is  the  same  frequency  as  the 
reference  frequency  FI.  If  the  divide-by-N  circuit  output 
frequency  is  initially  lower  than  the  reference  frequency  FI , 
the  oscillator  output  frequency  is  automatically  increased 
until  FI  and  FI*  are  identical.  When  this  occurs,  the 
oscillator  output  frequency  will  become  stable  and  remain  at 
that  frequency  until  the  number,  N,  changes.  When  this 
number  changes,  the  oscillator  frequency  is  again 
automatically  adjusted  so  that  FI  and  FI*  are  again 
identical. 

B.  PROPOSED  USES  OF  PROGRAMMABLE  FILTER 

A  wide  field  of  applications  exist  in  FH  systems  for  the 
programmable  GIC  filter.  The  outstanding  performance  of  the 
filter  (including  sensitivity  and  stability)  and  its  high 
speed  response  to  the  different  inputs  (due  to  the  use  of 
CMOS  integrated  circuits)  make  it  very  exceptional  in  this 
field.  The  first  proposed  use  is  indicated,  in  Fig.  (7.4). 
The  figure  illustrates  a  receiver  of  a  frequency  hopping 
demodulator.  The  received  frequency  hopped  signal  after 
heterodyned  by  the  RF  mixer  passes  through  the  GIC  filter  in 
a  BP  realization  (at  this  application  the  topology  control 
network  does  not  need  to  exist  since  the  BP  realization  is 


the  only  type  to  be  used).  The  frequency  shift  of  the  filter 
is  controlled  by  the  synchronized  pseudorandom  code 
generator.  This  code  generator  also  controls  Q  through  some 
interfaced  binary  logic  in  order  to  correct  the  amplitude 
reduction  which  appears  at  high  frequencies.  The  frequency 
shift  problems  can  be  easily  corrected  by  the  use  of 
composate  Op.  Amp.  Fig.  (7.5)  as  it  is  extensively  analyzed 
in  Refs.  [32],  [34],  and  [35]. 

The  programmable  filter  can  also  be  used  in  the  direct 
frequency  synthesizer  as  it  is  illustrated  at  Fig.  (7.2).  At 
each  frequency  mixer  output  exists  the  need  of  a  filter 
required  to  reject  the  unwanted  frequencies.  The  frequencies 
to  pass  are  not  always  the  same,  but  they  hop.  the  BP 
realization  of  the  filter  is  used  which  center  frequency  can 
be  controlled  accordingly.  The  filter  may  require  a  short 
time  for  the  signal  to  stabilize  after  it  is  selected.  The 
time  required  for  the  filter  to  •  be  stabilized  at  each  new 
frequency  may  ultimately  determine  the  maximum  hopping  rate 
of  the  direct  frequency  synthesizer. 


VIII.  CONCLUSION 


The  novel  design  described  here  has  resulted  in  a 
universal  programmable  filter  than  can  be  digitally 
controlled  to  realize  almost  any  practical  filter 
specifications.  This  is  done  through  the  use  of  CMOS 
switches  controlled  by  binary  codes  to  program  the  order  of 
the  filter,  the  filter  topology,  the  filter  center  frequency 
and  selectivity.  The  design  procedure  required  developing 
optimum  switching  arrangements  for  the  minimum  redundancy  in 
components  and  least  dependence  of  the  filtering  function  on 
switching  imperfections  such  as  switches  stray  capacitances 
and  non- zero  and  nonlinear  switch-on  resistance.  The 
sensitivities  of  Qp,  Wp  are  found  to  be  low  with  respect  to 
the  passive  and  active  elements  variations.  The  experimental 
result  show  close  agreement  between  theory  and  practice. 
Further,  these  results  indicate  that  these  realizations  are 
insensitive  to  temperature  and  power  supply  variations.  A 
wide  field  of  applications  exists  for  the  programmable  filter 
beside  the  one  discussed  in  Chapter  VII. 

(1)  Word  recognition  and  speech  synthesis; 

(2)  Music  applications; 

(3)  Signal  processing  in  communication; 

( 4 )  Adaptive  balancing . 


Further  investigation  is  needed  to  develop  a  prograiranable 
switched  capacitor  realization  that  can  allow  frequency 
scaling  by  changing  clock  frequency.  Work  can  be  also 
extended  for  developing  a  wide  bandwidth  programmable  filter 
using  the  composite  operational  amplifier  technique  proposed 
by  [39],  Such  implementation  would  lead  to  a  very  useful 
monolithic  device  at  moderate  cost. 

The  research  has  yielded  a  paper  that  was  presented  at 
the  19th  Annual  Asilomar  Conference  on  circuits,  systems  and 
computers,  Monterey,  California;  November,  1985  (Appendix  C). 
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APPENDIX  A 

t  JOB 

DIMENSION  ATI  SI  10  01  iPTISI  1001  lATZSt  1001 i PT2S(  1 00  I 
AT3S(  100)  •PT3SI  too ) »FAT 11  100) tFAT2( 100) trAr3{ 100) 

OI MENS I ON  ATI ( 100)«  PT II  too) »AT2(  100) •PT2I 1 00 ) • A 13 ( 1 00 ) t 
•  PT3I 1001 ,FPTl( 100) iFPT2| lOO) »FPT3l 100) i 

«  FTFI I  too) tFTF2ll00) «FTF3t 100 ) ,FP1I 100) ,FP2(  1 001 tFP3( t 00) 

COMPUEX  TlS|T2S|T33|YlStY2S»Y3S(YAStySS|y6StY7StYSStOS»KlStK2Sy 
«K3S,K ASfKSSfKdS 

COMPt.EX  TltT2»T3«Yt  »  Y2  •  Y3  »  YA*  Y5  *Y6»  Y7  ,  Y&  *  D  t  K 1  •  K2  t  K3  (KA  ,  K5  iKE 
DIMENSION  FRt 1001 tSTlSI 100) t^TtSl 100) t3T2S( 100) tHT2Sl 100) 

»tST3S( too) •HT3S( 100 ) fFBTt I  100) ,FBT2( 1001 «FBT3|  100) 

D I  MENS ION  BTl 1 100) , HTl I  100) »3T2( 100 ) »HT2( 100) (8731 100) 
ft«HT3f 100) ,FHT1| 100) tPHr2( 100) tFHT3( 100) 

COMPLEX  C1S,C2S |C3S*GlS«G2S»G3S«GAS*G5S»G6SyG7StGaS(00StXlS|X2S, 
«X3Sf XA3(XSS«  X6S 

COMPLEX  CltC2 »C3»Gt  «G2 • G3 * GA, GS»G6« G7 ,GS « OOy X 1 « X2 , X3 f  XA , XS«  X6,3 
OMEGA  s  0. 

DO  20  K  3  1,100 

omega  3  OMEGA  *  1 2S .920 . 

S  3  CMPLXI  0.0  ,  OMEGA  ) 

C 9999999999999999999  DATA  9999999999999999999999999999999999999999999 

R1  3  1600. 

Cl  3  2  0OE-O 
RQl  3  3200. 

R2  3  1600. 

C2  3  SOE-9 
Ra2  3  3200. 

R3  3  1600. 

C3  3  t2.5E-0 
RQ3  3  3200. 

RA  3  1600. 

CA  3  3.12E-9 
RQA  3  3200. 

A3  3  AA  3  A3S  3  AAS  3  lES 
A1  3  A2  3  AIS  3  A2S  3  lES 
J  3  A 
I  3  A 
L  3  A 
N  3  A 


9999999999999999SELECT  THE  TYPE  OF  THE  FOUPTH  FILTER  999999999999999 


DO  CASE 

L 

CASE 

MPF 

GIS 

s 

1  ./RA 

G2S 

s 

C  IS 

G3S 

s 

S9CA 

GAS 

s 

GIS 

G6S 

3 

GIS 

CSS 

3 

0.0 

GTS 

'  3 

C3S 

G8S 

S 

1  ./RQA 

CASE 

LPF 

GIS 

s 

l./RA 

G2S 

3 

S9CA 

G3S 

S 

G2SA  l./RQA 

GAS 

s 

GIS 

G63 

3 

0.0 

G5S 

3 

GIS 

GTS 

3 

0.0 

GBS 

xGlS 

CASE 

NOTCH 

GIS 

s 

l./RA 

G2S 

s 

GIS 

G3S 

s 

S9CA 

GAS 

s 

G  IS 

GSS 

s 

CIS 

G6S 

3 

0.0 

GTS 

3 

G3S 

I  ./R'34 


Gas  = 

1.-  NONE  OO 
C  BPP 

GIS  =  1./R4 
G2S  3  CIS 
63S  *  S«C4 
G4S  3  GIS 
GSS  3  0.0 
G6S  3  CIS 
G7S  3  1./Rq4 
Gas  3S»C4 
END  CASE 

c 

SELECT  THE  TYPE  OF  THE  THIRD  FILTER 

OO  CASE  N 


CASE 

HPF 


CASE 

LPF 


1./R3 

G1 

S«C3 

GI 

G1 

0.0 

G3 

t./RQ3 


1./R3 

S«C3 

G2  *  1./Rq3 


CASE 

NOTCH 


GI  3  I./R3 
G2  3  GI 
G3  3  S«C3 
04  3  OI 
G5  3  GI 
C6  3  0. 


C8  3  I./Rq3 
IF  NONE  OO 
C  8PF 

GI  31 ./R3 
G2  3  Gt 
G3  3  S#C3 
G4  3  GI 
G5  3  0.0 
C6  3  GI 
C7  3  i.yRQS 

oe  3S«C3 

ENO  CASE 

C  44444444444444P4SELCCT  THE  TYPE  OF  THE  SECOND  FILTER  «•**«« ft******** 
OO  CASE  I 
CASE 

C  MPF 

YIS  3  i./R2 
Y2S  3  YiS 
Y3S  3  S4C2 
Y4S  3  YIS 
V6S  a  YIS 
YSS  3  0.0 
Y7S  3  Y3S 
Yas  3  i./Rq2 


CASE 

LPF 


YIS  3  1./R2 


T3S  =  V2S 
Y4S  s  VIS 
Y6S  *  0.0 
Y5S  s  VIS 
Y7S  s  0.0 
vas  =YIS 

CASE 

NOTCH 

YIS  a  !.✓ 
V2S  a  YIS 
V3S  a  S*C 
VAS  a  YIS 
Y5S  a  riS 
Y6S  a  0.0 
Y7S  a  Y3S 
vas  a  1./ 
IF  NONE  OO 
3PF 

YIS  a  1 ./ 
V2S  a  YIS 
Y3S  a  S«C 
YAS  a  YIS 
YSS  a  0.0 
Y6S  a  YIS 
Y7S  a  !.✓ 
vas  aSftC2 
ENO  CASE 


Y2S«'  1./Ra2 
YIS 
0.0 
vis 
0.0 


1  ./R2 
VIS 
S«C2 
YIS 
YIS 
0.0 
Y3S 
1  ./RJ2 


1  ./R2 
YIS 
S«C2 
YIS 
0.0 
YIS 

1 ./RQ2 


SELECT  THE  TYPE  OF  THE  FIRST  FILTER 

OO  CASE  J 


CASE 

MPF 


CASE 

LPF 


1  ./Rl 
Yl 

S«C1 

VI 

Yl 

0.0 

V3 

i./RQl 


l./Rl 

S«Cl 

V2  *  l./RQl 


CASE 

NOTCH 

VI  a  I . /R  t 
V2  a  Yl 
Y3  a  SACl 
VA  a  Yl 
Y5  a  Yl 
Y6  a  O. 

YT  a  y3 
Y3  a  l./RQl 
IF  NONE  OO 
aPF 

Yl  a  l./Rt 
Y2  a  Yl 
Y3  a  S«C1 
YA  a  VI 
Y5  a  0,0 
Y6  a  Yl 
Y7  a  l./RQl 
Va  aSPCl 
ENO  CASE 


c 

FRIK)  s  QMcGA/6«28 
S  -  CMPLXI O.Of OMEGA) 

A13  =  ris«|l.»  1./A2S) 

K2S  =  I y2S  ♦  Y5S  ♦  Y6SI 
K3S  =  JY4S  ♦  Y7S  ♦  YSS) 

KAS  3  YAS  •»  Y8S 
KSS  3  YIS  A  Y3S 
K6S  3  YSS  A  Y6S 
C 

K1  3  Y1«U.A  t./A2) 

K2  3  tY2  A  YS  A  y6) 

K3  3  IYA  A  Y7  A  Y8I 
KA  3  YA  A  Y8 
K5  3  Yi  «  Y3 
K6  3  YS  A  Y6 
C 

Xts  3  GISAM-A  l./AASl 
X2S  3  tG2S  A  G5S  a  G6S) 

X3S  3  (GAS  A  G/S  A  G8S) 

XA3  3  gas  a  G8S 
XSS  3  GIS  A  G3S 
X6S  3  GSS  A  G6S 
C 

XI  3  Gl«(t.A  l./AA) 

X2  3  (G2  A  GS  A  G6) 

X3  3  fGA  A  G7  A  G8) 

XA  3  GA  A  GS 
X5  3  GI  A  G3 
X6  3  GS  A  G6 
C 

c 

OS  3(KtS«K2S«<3S)/AlS  A  YtS*YAS«K6S  AY3S3K2SAK3S/ 1 A  1 SPA2S) a 

*  Y2S«Y3S«YAS/A2Sa( Y2S/YtS)«KlS«Y3S«(Y7SAY8St  A  Y3S«K6S«K3S/A 2S 
T1S3(YIS«YAS#YSS  aY3S#V5S*(Y7S-Y8S)  aY3S*Y7SA( Y2Sa ySS ) a ( Y7 5« 

*  (Y1SaY3SI«K2S)/AIS)/OS 

T2S3(YIS«YSS«KAS  -YIS*Y7S#Y6S  A Y3S#Y5S#KAS/A2S  -Y3S*Y6S#Y7S/A23 
A  AY5SAY7SAKSS/AIS  A  Y2SAY3SAY7S  ) /OS 

T3S3 (KlS«Y7S«K2S/AtS  AKtS9YAS«YSS  a Y3SAY7SAK2S/A tS*A2S  AY3S9YAS 
«  9Y5S/A2S  A  I YtS/Y2S}«KlS*Y3S«Y7S  a  K6S« Y3S9Y7S/A2S ) /OS 
C 

O  S(K1«K2«K3)/At  a  Y19YA9K6  aY39I<29K3/!  A19A2)  A  Y29Y39YA/A2  A 

*  <Y2/VI >«K1*Y3*I Y7AY8)  a  Y39K69K3/A2 

Tl 3( Y1«YA«YS  A Y3AYS*( Y7-V8I  AY39 Y7»l Y2aY6| A C  Y7A ( Y I AY3 ) AK2 ) /A 1 J /O 
T23(YI*Y5*KA  -V19Y79Y6  AY3#Y59XA/A2  -Y39Y69Y7/A2  AY59Y79K5/A2  A 

*  Y29Y39Y7  I/O 

T33(KI*Y7AK2/A1  aK 1*YA«YS  AY39Y70K2/AI9A2  aY3»YA9Y5/A2  a(Y2/YI)9 

*  K19Y3AY7  aKG4Y39Y7/A2)/0 

C 

c 

OOS  3( X1S«X2S9X3S}/A3S  a  GtS*GAS*X6S  AG3SftX2S«X3S/( A3S9AAS) a 

*  G2S*G3S«GAS/AAS A(G2S/GtS) *XlS«G3S«tG7SAGasl  a  C3S9X6S9X3S/AAS 
Cl S34G1S«GAS«GSS  AG3S9GSS9IG7S-C8S)  aG3S»G7S«(G2SaG6S) a( G7S9 

»  (G1SaG3S)«X2SI/A3S)/OOS 

C2S3(G1S«G5S9XAS  -G1S9G7S9G6S  AG3S9GSS9XAS/ AAS  >G3S*G6S9G7S/AAS 
«  AG5S9G7S9XSS/A3S  A  C2S9G3S9G7S  )/OOS 

C3S3(xiS«G7S«X2S/A3S  aX1S«GAS«G5S  aG3S»G7S«X2S/A3S*AAS  aG3S*GAS 

*  ftGSS/AAS  A  I G1S/G2S}«X1S4G3S«G7S  a  X6S9G3S4G7S/AAS )/00S 
C 

OOsf X1*X2*X3)/A3  a  G19GA9X6  aG3«X2«X3/< A3*AA| a  G29G39GA/AA  A 
«  (G2/G1 )«X1 «G3«|G7aG8)  a  G3«X6aX3/AA 

Cl3( Gl«GA«GSAG3«GS«(G7-Ga)AG3*G7«(G2AC6> A( C7*|G1aG31«X2)/A3)/D0 
C23( G1«G5*XA-G1«C7«G6  aG3«GS*XA/AA  -G3«C6*G7/AA  aGS«G7«XS/AA  a 

*  G2«G3«67  1/00 

C33(X1«G7«X2/A3  aX1*GA«GS  a'.3«G7«X2/ A39AA  aG3«GA«GS/AA  a(G2/G1|« 

*  X1«G3«G7  AX8«G3«G7/AA 1/00 


ATlSfK)  3  20.9A1.0G10ICABSIT1S)) 
AT2S(K)  3  20  .aALOGlOICABSI  T2Sn 
AT3S(K)  3  20.«Ai.0G10ICABSf  T3Sn 


ARTIS  s  REAL(TIS) 

ART2S  a  REAL|T2S> 

ART3S  a  REALIT3S) 

AITIS  a  A(MA3(  T1  S) 

AIT2S  a  AIMAG<T2S} 

AIT3S  a  AIMAGtT3Sl 

PT1S<<)  a  ArAN2<AXTtS  t  ARTtS}»  57.325 
PT2S(K)  a  ATAN2IA1T2S  »  ART2S)tt  57.325 
PT3SIK)  a  ATAN2«A1T3S  ,  ART3Sl«  57.325 

ATlIKI  a  20.«A|.aGI0lCAeS(Tl  II 
AT2tKl  a  20.«AUOG10(CAaS(T2}| 

AT3fK|  a  20.«ALaG10<CAaS(T3l I 

ARTl  a  REALl TX 1 
ART2  a  REAI_(  T2  I 
ART3  a  REAI.IT3I 

AITI  a  AX»AG(TX} 

AIT2  a  AIMAG1T2I 
AXT3  a  AIHAG4r3l 

PTHKI  a  ATAN2fAXTl  ,  ART!  |»  57.325 
PT2(K)  a  ATAN21AXT2  »  ART21«  57.325 
PT3rK|  a  ATAN2IAXTJ  ,  ART3)«  57.325 

FATIIKI  a  ATX4KI  ♦  ATlSXKl 

FAT2IKI  a  AT2(KI  *  AT2SIKI 

FAT31K}  a  Ar3<X|  «■  AT3StK) 

FPT14KI  a  PTlfKl  ♦  PTlStKt 

FPT2t<l  =  PT2IKI  ♦  PT2SIKI 

FPT3IKJ  a  PT3fK|  ♦  PT3S4K) 

ariSIKI  3  20.PALOGXOiCASSlC15l I 
5T2S4K}  a  20 .PAUOC 10 ( CASSf C2S} I 
aT3S4K|  a  20.PAL0C104CABS4C3SI ) 

ARTIS  a  REAL4CXSI 
ART2S  a  REAI.4C2S} 

ART3S  a  REAL4C3S} 

AITIS  a  AXMAGtClSl 
AIT2S  a  AXMAGIC2S} 

AXT3S  a  AIMAG4C3S1 

MTIS4K)  a  ATAN2IAXT1S  »  ARTXSIP  57.325 
HT2S4KI  a  ATAN2tAIT2S  ,  ART2SI4  57.325 
HT3S4K)  a  ATAN24A1T3S  ,  ART3Si«  57.325 

3T14K)  a  20.«AU0GI04CAaSICl}| 

8T24KI  a  20.«ALaGI04CABSIC2}} 
aT3IKl  a  20.*AI_aG104CA8S(C3|} 

ART!  a  REALlCII 
ART2  a  REACf  C2 I 
ART3  a  REAI.XC3I 

AITI  a  AXMAG4C11 
AXT2  a  AiMAa4C2l 
AIT3  a  AIMAG4C31 

HTtIKl  a  ATAN24AlTt  i  ARTl |*  57.325 
HT2fK|  a  ATAN24AXT2  •  ART2|«  57.325 
HT3fK}  a  ATAN2fAIT3  (  ART3|«  57.325 

FBTltKI  a  STIIKI  ♦  BTlStKl 

F8T24K)  a  8T2IKI  ♦  BT2SfKl 

FeT3|K}  a  BT34K}  ♦  BT3S4Kt 


tPKOU 


EHROU 


%JOS 


APPENDIX  B 


j 

0 

j 


■J 

3 


X 


DIMENSION  ATI  SI  100}  fPTISl 100} tAT2Si  100) tPTZSf  100) 

«»  AT3S(  100)  (PT3S  (  t  OO)  tFATK  tOO)  fFAT2t  100)  (FATB  (100) 

OI  MENS  ION  ATI  1  too)  .PTK  100)  fATaC  100)  »PT2I  100)  tAT3(  100)  t 
«  PT3( 100) tFPTlI  too)  (FPral too) »FPT3(  100) t 

»  FTFl  1100)  »FTF2(  100  )  »FTF3t  100)»FP1(  100)  t F=P2 1 1  00 )  * FP3 1  1  00 ) 

COMPLEX  TlStTaSf T3S  tV IS  * Y2S «V3S » YAS »ySS , VoSt V7S » V3S ,OStKlS»K2S* 
«K3StXAS»XSS|K6S 

COMPLEX  T1,T2»T3«Y1 « Y2« Y3 • YA* YS»Y6» Y7 » Ya »0^ K1 • K2 t K3 (KA, KS « K6 
DIMENSION  FRI  lOO)  *3 T ISl 1 00) iHT ISl 1001  »aT2S( 1 00 1 1 HT2S ( 100 ) 

AtBTSS  (  too)  *HT3S(  100  )  fFBTK  100)  tFBTai  100)  *FBT3(  100) 

DIMENSION  OTl (1 00),HT1| 100) tBT2l 1 OO ) , HT2 I  100) (BT3( 100) 

«tHT3(  100) (FHTII  lOO)  (FHT2( 100)  »FHT3(  100) 

COMPLEX  ClS(C2S«C3S«ClS(G2S»G3S«GAS»6SSyG6S»G7St63S,OOS«  X1S|X2S« 
»X3S ( X AS  t XSS  t  X6S 

COMPLEX  C1«C2«C3«G1 iGatCa * GA« G5«G6t 07 (GO *00 i X 1 •X2 *X3 » XA , X5« X6 tS 
OMEGA  3  0. 

OO  20  K  3  1,100 

OMEGA  =  OMEGA  *  12S.A20. 

S  3  CMPLXI  0.0  ,  OMEGA  ) 

DATA  «««««««« AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

PI  3  1600. 

Cl  3  200E-9 
RQl  3  3200. 

R2  3  1600. 

C2  3  SQE-9 
RQ2  3  3200. 

R3  3  1600. 

C3  3  12.SE-9 
RQ3  3  3200. 

RA  3  1600. 

CA  3  3.12E-9 
RQA  3  3200. 

Wl  3  wis  3  "SET  VALUE  " 

A3  3  AA  3  A3S  3  AAS  3  Vl/S 

At  3  A2  3  AIS  3  A2S  3  WIS/S 

J  3  A 
X  3  A 
L  3  A 
N  3  A 


c  ft««33«««3«3*««««SELECT  The  type  of  the  FOURTH  FILTER 


3  l./RA 

3  GIS 
3  S«CA 
3  CIS 
3  GIS 
3  0.0 
3  G3S 
3  1 ./RQA 


3  l./RA 

3  S«CA 

3  G2SA  l./RQA 
3  GIS 
3  0 .0 
3  CIS 
3  0.0 
3GtS 


3  l./RA 
3  CIS 
3  S*CA 
3  CIS 
3  CIS 
3  0.0 


OO  CASE  L 
CASE 
HPF 

GIS 

G2S 

G3S 

GAS 

G6S 

GSS 

GTS 

Gas 

CASE 

LPF 

GIS 

G2S 

G3S 

GAS 

G6S 

GSS 

GTS 

Gas 

CASE 

NOTCH 

GIS 

G2S 

63S 

GAS 

GSS 

G6S 


185 


.r 


68S  3  t  ./'R'34 
NONE  OO 

C  BPF 

GIS  3  1./R4 
G2S  3  GIS 
G3S  3  S«C4 
G4S  3  GIS 
GSS  3  0.0 
G6S  3  GtS 
G7S  3  t./Rq4 
GSS  3S«C4 
END  CASE 
C 

C4949944944444444  SELECT  THE  TYPE  OP  THE  THIRD  FILTER 
OO  CASE  N 
CASE 

C  HPF 

GI  3  1./R3 
G2  3  GI 
G3  3  S«C3 
G4  3  GI 
G6  3  GI 
G5  3  0.0 
G7  3  G3 
G8  3  t./RQ3 

CASE 

C  LPF 

Gt  3  I./R3 

G2  3  S«C3 

G3  3  G2  •  1./Rq3 

G4  3  GI 

G6  3  0.0 

GS  3  GI 

G7  3  0.0 

G8  3  GI 

CASE 

C  NOTCH 

GI  3  I./R3 
G2  3  GI 
G3  3  S*C3 
G4  3  GI 
GS  3  Gt 
C6  3  0. 

G7  3  G3 
G8  3  I./Rq3 
IF  NONE  OO 
C  8PF 

GI  3  I./R3 
G2  3  GI 
G3  3  S«C3 
G4  3  GI 
GS  3  0.0 
G6  3  GI 
C7  3  1./Rq3 
G8  3S*C3 
END  CASE 

C  4««4ft«««»««««««9SELCCT  THE  TYPE  OF  THE  SECOND  FILTER 
OO  CASE  I 
CASE 

C  HPF 

YIS  3  I./R2 
Y2S  3  YIS 
YJS  3  S4C2 
Y4S  3  YIS 
Y6S  3  YIS 
YSS  3  0.0 
Y7S  3  YSS 
Yes  3  1./RQ2 

CASE 

C  LPF 

VIS  3  1./R2 


I  L>  Ai  tjUV 


Y3S  s  V2S>  l./R<32 

Y*S  =  YIS 

Y6S  =  0.0 

Y5S  =  YIS 

Y7S  s  0.0 

ras  srts 


CASE 

NOTCH 

YIS  a  1./ 
Y2S  a  YIS 
Y3S  a  S*C 
Y4S  a  YIS 
YSS  a  YIS 
Y6S  a  0.0 
Y7S  a  Y3S 
YSS  a  1  .y 
IF  NONE  OO 
BPF 

YIS  a  1 ./ 
Y2S  a  YIS 
Y3S  a  s«: 
Y*S  a  YIS 
YSS  a  0,0 
Y6S  a  YIS 
Y7S  a  l.y 
YSS  aSftC2 
ENO  CASE 


t  ./R2 
YIS 
S«C2 
YIS 
YIS 
0.0 
Y3S 
1  .yR22 


1  ./R2 
YIS 
S«C2 
YIS 
0.0 
YIS 

1 ./RQ2 


SELECT  THE  TYPE  OF  THE  FIRST  FILTER 

OO  CASE  J 
CASE 

C  HPF 


CASE 

LPF 


CASE 

NOTCH 


Yl 

a  l./Rl 

Y2 

a  Yl 

Y3 

a  S#Cl 

YA 

a  Yl 

Y6 

a  Yl 

Y5 

a  0.0 

Y7 

a  Y3 

Y8 

a  l./RQl 

Yl 

a  l./Rl 

Y2 

a  S«Cl 

Y3 

a  Y2  ♦  t./ROl 

YA 

a  Yl 

Y6 

a  0.0 

Y5 

a  Yl 

Y7 

a  0.0 

Y8 

a  Yl 

1 

Yl 

a  l./Rl 

Y2 

a  Yl 

Y3 

a  S«C1 

YA 

a  Yl 

Y5 

a  Yl 

V6 

a  0. 

Y7 

a  Y3 

Y3 

a  l./RQl 

IF  NONE  OO 
SPF 

Yl  a  t./Rl 
Y2  a  Yl 
Y3  a  S*Cl 
YA  a  Yl 
Y5  a  0.0 
Y6  a  Yl 
Y7  a  l./RQi 
Y8  aS«Ct 
ENO  CASE 


UU  V  W  «J  uu 


c  «*««*«*««««»«*««««**««  «««»»»«»« 

c 

FniK)  =  □MEGA/6.2e 
S  =  CMPLXt a*0»aMEGAl 
AIS  S  1./A2S) 

K2S  s  ( Y2S  »  VSS  *  V6S) 

K3S  =  (Y4S  ♦  r7s  ♦  vas) 

KAS  a  Y4S  ♦  Y8S 
KSS  a  YIS  ♦  Y3S 
K6S  a  VSS  »  Y6S 

c 

K1  a  Vl«(1.4  1./A2) 

K2  a  |V2  ♦  VS  ♦  Y6l 
K3  a  fVA  ♦  Y7  ♦  VSl 
K4  a  Y4  ♦  V8 
K5  a  VI  ♦  V3 
K6  a  VS  ♦  V6 
C 

XIS  3  Gis*(l.»  1./A4SI 
X2S  a  iG2S  *■  G5S  ♦  G6St 
X33  a  IG4S  *  GTS  *  GaS) 

X43  a  G4S  *  G8S 
XSS  a  CIS  »  G3S 
X6S  a  GSS  ♦  G6S 

c 

XI  a  Gl«(t.*  t./A4) 

X2  a  (G2  ♦  GS  ♦  G6| 

X3  a  (G4  ♦  G7  »  Ga) 

X4  a  G4  »  G8 
XS  a  G1  «■  G3 
X6  a  GS  *  G6 


OS  a<KlS«K2S«<3SI/AlS  *■  VtS*Y4S»X6S  *-V3SCK2S»K3S/ f  A  I  S»A2S)  * 

•  Y2S*Y3S*V4S/A2SA<  Y2S^YlSI*iKlS*Y3S4lY7S>Y8S»  ♦  Y3S#K6S4«C3S/A  2S 
TlSa4YlS*Y4S«YSS  ♦ VSSaYSS* I Y7S- YSS )  ♦Y3S4Y7S4f Y2S*Y6S} ♦( Y7S4 

*  (YISAY3SI*K2SI^AlSi/OS 

T2Ss4YIS«YSS«X4S  -YIS*Y7S«Y6S  ♦  Y3S*Y5S*IC4S/A2S  -Y3S4Y6S*Y7S/A2S 

*  ♦YSS«Y7S*KSS/A1S  ♦  Y2S*Y3S*Y7S  >/OS 

T3Sa  (K1S#Y7S«K2S/A1S  ♦■KIS»Y4S*Y5S  ♦Y3S4Y7S4K2S/A  ISO A2S  ♦Y3SOY4S 

*  OYSS/A2S  ♦  I  YtS/Y2SJ«ltIS*V3S*Y7S  ♦  K6SOY3SOY7S/A2S ) /OS 

O  a|Kl*K2*K3|/Al  ♦  YI«V4ftK6  ♦  Y3««K2*K3/t  AIOA2 )  ♦  Y20Y30Y4/A2  ♦ 

»  IY2/YI >*KIOY3*l Y7»Y8>  ♦  Y30K60K3/A2 

TIa(Yt*Y40Y5  ♦ Y34YS* « V7-Va I  aYSO Y7»|  Y2^Y6) ♦! Y7*| V I ♦YS ) OK2 ) /A  I ) /O 
T23<VI«Y5»K4  -Y10Y70V6  ♦  V30Y  50X4/A  2  -Y30Y60Y7/A2  ♦YS0V70K5/A2 

♦  Y20Y30Y7  )/0 

T3a« KI0V70K2/A 1  AK I0Y40Y5  ♦Y30Y70K2/AI0A2  ♦Y30Y40Y5/A2  ♦IY2/Y1)* 

•  Kl#Y30Y7  ♦KSOV3OY7/A2I/0 


□OS  af X1SOX2SOX3S1/A3S  4  GIS«G4S«X6S  4G3SftX2S«X3S/| A3S4A4S} 4 
G2S«G3S«G4S/A4S4(G2S/GtSl*XlSOG3S*(G7S4C8S)  4  G3S*X6S«X3S/A4S 
Cl Sa(GtS«G4S*GSS  4G3S9GSSO |G7S-GaS )  4G3SoG7S«|G2S4G6S1 4| C7sa 
IG1S4G3S}*X2SI/A3S}/OOS 

C23a(GIS*GSS9X4S  'G I SOG7S9G6 S  4G3SOGSSOX4S/A4S  •G3S9G6S9G7S/A4S 
4G5S9G7S9XSS/A3S  4  G2SOG3SOG7S  )/DOS 
C3Sa|xiS«G7S«X2S/A3S  4XtS«G4S«G5S  4G3SftG7S*X2S/A3S«A4S  4G3SOG4S 
9G3S/A4S  4  I G1S/G2S10X1S0G359G7S  4  X 6S9G3S9G7S/A4S )/OOS 

DOaf XI 9X29X3 1/A3  4  G19G49X6  4G39X29X3/< A39A4| 4  G29G39G4/A4  4 
tG2/Gl )9Xt 9G39tG74Gai  4  G39X69X3/A4 
Ct3( G19G49GS4GJ9G59<G7-Ga)4G39G79lG24G6)4t  G79l G14G3) 9X21 /A3 i /OO 
C2a| 619GS9X4-Gt9G79G6  4G39G59X4/A4  -G39G69G7/A4  4GS9G79X5/A4  4 
G29G39G7  | /OO 

C3a|X19G79X2/A3  4X19G49GS  4G39G79X2/ A39A4  4G39G49GS/A4  4CG2/G1)9 
X19G39G7  4X69G39G7/A4I/00 


ATISfX}  a  20.9A1.0G10f  CABSI  TIS)  ) 
AT2SIK)  a  20  .OALOG 104 CABSI  72511 
AT3SfK)  a  20.«At.(X:iO(CABS{  T3Sn 


c 

ARTIS  3  REALiriSl 
ART2S  =  REAHTaS) 

ART3S  3  AEAL(T3S| 

C 

AI  TIS  3  Af  MAGI  T1  SI 
AIT2S  3  AIMAG<T2S) 

AIT3S  3  AIMAGIT3SI 
C 

PTISIKJ  3  ATANEIAZTIS  »  ART>S)»  S7.32S 
PT2SIK>  3  ATAN2IAIT2S  »  ART2S)#  S7,32S 
PT3SIK)  3  ATAN21A1T3S  •  ART3S)«  S7.32S 

^  ATllK)  s  20.»Al.aGt0<CAaSf  Tt  II 

AT2IKJ  3  20.«ALOGta(CA8S(T2n 
AT3IK1  3  20.9AI.OG104CA8StT3)  I 
C 

ARTl  3  REALITl) 

ART2  3  REAi_|T2» 

ART3  3  REAI_<T3| 

C 

AITl  3  AI mag  I  Til 
AtT2  3  A1MAG(T2) 

AIT3  3  A1HAGIT3) 

C 

PTllKl  3  ATAN2|AZTt  •  ARTli*  57,325 
PT2tK)  3  ATAN21A1T2  ,  ART2)«  57.325 
PT3tlCI  3  ATAN2(AITJ  ,  ART3)«  57.325 
C 

FATt(K)  3  ATllKl  ♦  ATtSIKl 
FAT2(K|  3  ATalKI  ♦  AT2StK) 

FAT3(K)  3  AT3<KI  «•  AT3SIIH 
C 

FPTIIK)  3  PTZiKl  ♦  PTIStKI 
FPT2IKI  3  PT2fKl  ♦  PT2SIKI 
FPTaiKl  3  PT3<Ki  ♦  PTSSIKI 

STISIK)  3  20.«AI.QGt0(CABSICZSl  I 
aT2SIK)  3  20.9AI.0Gt0(CAaSIC2S)  ) 
aT3S(K)  3  20«3ALQC10(CA8SIC3Sn 

c 

ARTIS  3  REALICIS) 

ART2S  3  REAL(C2S1 
ART3S  3  REAL(C3S1 
C 

AITIS  3  AIMAGiCtSI 
Alias  3  AIMACfC2S} 

A1T3S  3  AIMAG(C3S1 
C 

HTISIK)  3  ATAN2fAITtS  $  ARTtS)«  57.325 
HT2SIKI  3  ATAN2(A(T2S  •  ART2S|«  57.325 
HT3S(K|  3  ATAN2IA1T3S  •  ART3S)*  57.325 
C 

3TIIKI  3  20.«AI.0Gt0lCABSlCl  >) 

BTalK)  3  20.«Al.aG10<CABSIC2l) 
aT3(Kl  3  20.«AI.OGialCABSIC3}  J 

c 

ART!  3  REAUtCl) 

ART2  3  REAI.IC2) 

ART3  3  REALIC3) 

C 

ATTt  3  AIMAGIClt 
Aira  3  aimag(C2) 

A1T3  3  A1NAG(C3) 

c 

HTtIKl  3  ATAN2IAIT1  t  ARTl I*  57. 3ZS 
HT2<ICl  3  ATAN21ATT2  •  ART2)*  57.325 
MT3IK)  a  ATAN2fAIT3  »  ART3I*  57.325 
C 

FBTIIKI  a  3TIIKI  ♦  BTISIKI 

FBTaiK)  a  BTZIK)  ♦  BT2SIK1 

FBT3IKI  a  BT3fK|  ♦  BT3SIK} 
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FMTl tK» 
FMT2 IK) 
FMT3IK ) 

FTFl IK) 
FTF2IK) 
FTF3I<) 


MTliKl  ♦  V 
MT2IK)  ♦  V 
HT3IK)  ♦  » 

ATI  IK)  ♦ 
FAT2IK)  ♦ 
FAT3IK)  ♦ 


FPt(K)  =  FHTlfK) 
FP2IK)  a  FHT2<K) 
FP3IK)  s  FMT3IK) 


HTlSf K) 
HT2SI K) 
HT3SIK) 

STlfK)  » 
►  F3T2IK) 
FBT31K) 

FPTICK) 

FPT2IK) 

FPT3IK) 


eTlSlK)«  ATISIK) 


♦ENTRY 


WRITE! 6,66)  FR(K }«FATl|K),F3Tt |K ) 

*  »FTF1|K) yHTt|K),FHri (K) 

»  FHT3IK) 

CONTINUE 

FQRMATI6I 1X(F9.3) ) 

STOP 

END 
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ABSTRACT 

In  tnia  contribution  a  general  purpose 
digitally  controlled  analog  filter  Is  presented. 
The  novel  design  la  a  cascade  of  second  order 
sections  that  are  Individually  programmed  to 
achieve  any  filtering  topologies.  TwO‘~binary 
words  are  used  to  control  the  pole  frequency  up 
and  selectivity  Qp  of  each  section  Independently. 
Each  second  order  section  is  a  Generalized* 
Immlttance  converter  (GIC)  biquads  which  are 
known  for  their  high  stability  and  low  active 
and  passive  sensitivity.  CMOS  switches  are  used 
to  electronically  relocate  the  minimum  number  of 
passive  elements  to  achieve  function 
programmabt  11  cy  .  Switches  are  also  used  to 
select  the  number  of  cascaded  sections  to 
realize  higher  order  transfer  functions. 

INTRODUCTION 

The  availability  of  an  analog  filter  with 
digitally  controlled  "programmable"  coefficients 
has  been  the  goal  of  many  researchers  due  to  Its 
several  attractions.  One  possibility  of  a  com* 
pact,  versatile  analog  filter  under  remote 
control  opens  up  many  novel  and  Independent 
application  areas.  Also,  when  a  programmable 
filter  la  combined  with  a  permanent  reference 
memory  which  is  user*programniable,  this  would 
form  an  economical  and  versatile  device  for 
dedicated  stand*alone  applications.  The  need 
for  such  a  device  was  motivated  by  advancement 
in  film  and  semiconductor  technologies  as  well 
as  the  continuous  upgrading  of  systems  spectflca* 
tions  to  tsk#  advantege  of  the  available 
technologies  to  tne  limits. 

Linear  analog  filtering  finds  many  applies*' 
tions,  such  as  speech  processing  (recognition  or 
synthesis!,  geology,  Instrumentation,  coomuntes* 
tions,  process  control,  adaptive  balancing,  etc. 
There  has  been  much  emphasis  on  performing  the 
filtering  function  digitally,  largely  because  of 
ins  ease  of  varying  and  optimizing  the  transfer 
function.  However,  for  many  reasons,  such  as 
cost,  slzs,  signal  processing  complexity,  and 
bandwidth.  It  would  be  desirable  to  perform  the 
filter  function  with  linear  components  yet 
retain  tne  flexabtlity  of  varying  the  filter 
parameters  digitally. 

Recently,  the  advantages  of  combining 
linear  oomponents  (operational  amplifiera  (OAS), 


resistors  and  capacitors)  and  nonlinear  elements 
(switches)  have  been  demonstrated  using  switched 
capacitor  techniques  In  this  contribu¬ 

tion.  we  are  presenting  the  results  of  realizing 
a  continuous  active  device  using  linear  elements 
and  switches  controlled  by  digital  signals  to 
achl.eve  a  fully  programmab^t  filter  [4]. 
Several  programming  features  of  the  proposed 
filter  are  reported.  The  first  feature  is  the 
ability  of  the  network  to  realize  the  most 
common  filtering  functions  (function  programma¬ 
bility)  namely:  Low  Pass  (LP).  Band  Pass  (BP), 
High  Pass  (HP).  All  Pass  (AP)  and  Notch(N) 
functions,  using  the  minimal  set  of  elements. 
Tne  second  feature  is  the  ability  of  the  network 
to  program  (Independently)  the  key  parameters  of 
the  filtering  function  cnosen  (parameter  program¬ 
mability)  namely:  the  pole  resonent  frequency 
(up)  and  selectivity  (Qp).  Finally  the  ability 
to  program  the  network  to  cascade  several  sec* 
tions  to  achieve  higher  order  filter.  All  of 
the  above  programmability  features  are  performed 
independently  to  realize  a  universal  filtering 
network . 

2.  DESIGN  ANALYSIS  OF  THE  PROPOSED  FILTERS 

The  basic  active  network  considered  here  as 
the  heart  of  the  programmable  filter  is  the 
second  order  Generalized  Immlttance  Converter 
(OlC)  structure  T^j,  Fig.  1,  wnose  superior 
performance  was  established  In  tne  literature 
[6].  The  general  transfer  function  realized 
T(5)  is  given  by: 

T(s)  -  M(s)/D(s)»(ap*a<s*as2)/(bQ*b^s*b2S^)  (1) 
The  ClC  transfer  functions  of  Fig.  1  assuming 
non-ideal  OAs  are  given  by 


*1  -'r’l 

'l'*’ 

t,... 

0  ,, 

•T, 

’J 

Stability  and  sensitivity  analysis: 


An  important  criterion  of  a  realization  la  its 
sensitivity  to  slsment  variations.  The  ClC 
sensitivity  analysis  has  shown  to  be  as  good  or 
bettsr  than  all  competitive  second  order  net¬ 
works  (6}.  While  the  GIC  stability  can  easily 
bf  demonstrated  since  in  all  the  transfer 
functions  (2),  the  coefficients  of  D(s)  art  seen 
to  remain  positive  for  any  OA  mismatch.  This  is 
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dj«  to  ttit  absene*  of  nagatlva  taraia  In  D(a). 
Thtrofore  tna  zaros  of  D(a)  will  raaain  tn  tna 
laft-nair  s-plana  and  low  fraquoney  unitabla 
■oOas  cannot  arlsa  during  aetlaatlon. 

Function  prograaaaabllltyi 

Tna  objaetlva  of  thla  raaaareh  waa  to 
davalop  a  daalea  that  la  capabla  of  raallzlng 
the  following  tranafar  functional  LF  wnara  T(a) 
•  K/0(si.  BP  wnara  T(a)  •  KS/D(a).  HP  wnara  T(a) 
-  Ks2/D(a).  aP  wnara  Ttal-la^-alwp/Qpl^Wp^l/Ola) 
and  N  wnara  T(a)>(a2<ux^)/0( a) ,  By  optlalzlng 
tna  daaign  of  tha  flltar.  It  waa  found  tnat  all 
of  tna  abova  functlona  can  ba  raallxad  by  tna 
aacond  ordar  QIC  aactlon  uatng  four  raalatora, 
two  capacltora  and  two  OAa  aa  anown  In  Table  1. 
Tnaaa  paaalva  alcaanta  arc  eonnactad  to  dlffar- 
ant  nodaa  to  acniava  tna  varloua  raallzatlona. 
A  aat  of  CMOS  bllataral  awltcnaa  eontrollad  by 
digital  binary  word,  ara  uaad  to  ralocata  tna 
aana  alaaanta  In  dlffarant  waya  to  acniava  tha 
daalrad  filtering  functlona  according  to  Fig.  2. 
Tha  truth  table  of  the  awltchaa  control  logic. la 
anown  In  Table  2,  wnila  Fig.  3  llluatrataa  tha 
corraaponding  alnlalzed  CHOS  logic  circuit  uaad 
for  paaalva  alaaenta  relocation. 

Paraaatar  prograaviabllltYi 

Uhl  la  four  of  tna  raalatora  ara  aqual  and  of 
value  R  each,  the  fifth  raalator  la  the  Qp 
datarainlng  raalator  and  of  valua  Rq*RQp.  Tna 
two  capacltora  are  aqual  and  of  valua  e«i/(upR) 
aacn.  Two  equal  banka  of  capacltora  ara  uaad  to 
control'wp.  Caen  bank  oontalna  n  binary 
waigntad  capacitora  connaetad  In  aarlas  with  n 
CH0;S  awltchaa  aa  anown  In  Pig.  A.  Ualng  n  bit 
binary  word  to  control  tna  awltcnaa,  2" 
different  valuaa  of  o  can  ba  obtained  that  eor- 
rasponda  to  2"  dlffarant  valuaa  of  wp.  Ualng  a 
alallar  tacnnlqua,  tna  valua  of  Rq  can  ba  con* 
trolled  by  an  a  bit  digital  word  tnat  ylalda'2* 
different  valuaa  of  Qp,  aa  llluatratad  In  Fig.  S. 
Thua,  full  indapandant  control  of  the  pole  pair 
■In  and  Qp  ara  aenievad  by  prograanlng  the 
digital  worda  controlling  the  awltcnaa  to  obtain 
tna  corraaponding  o  and  Rq.  The  coaplata  aacond 
order  prograaaabla  flltar  la  anown  In  Flg.£.|, 
wnara  tna  funetlon  programabll i ty  aa  wall  aa 
the  paraaatar  prograaaabl 1 I ty  ara  daaonatratad. 

Higher  ordar  proaraaaabllltr: 

Actlva  flltera  daaign  prooadura  can  bn 
claaalfind  aa  direct  or  eaaeade.  In  dlraot 
ayntneala  proeaduraa  the  tranafar  funetlon  la 
realized  aa  a  aingla  aootlon  CT].  In  eaaeada 
ayntneala  proeaduraa  a  high  ordar  tranafar 
function  la  ezpreaaad  aa- a  product  of  flrat  and 
aacond  ordar  tranafar  functlona  and  each  of 
tneae  la  realized  Independently.  Tna  overall 
network  la  obtained  by  eaaeading  tna  Individual 
aaetlona.  The  eaaeade  aathod  of  aynthoata 
offera  two  practical  advantagaa  (a)  alapla 
network  tuning  (b)  a  few  nuabor  of  unlvaraal 
aaetlona  can  ba  dealgnad  which  can  raallta  a 
aultltuda  of  network  apaolfloatlona. 


The  aacond  ordar  QIC  network  atructura 
lands  itself  to  tna  eaaeada  aynthasls  procedure 
since  It  does  not  raqulra  additional  Isolating 
aapllflara.  Flg.6.b  shows  a  block  dlagraa  of  a 
prograaaabla  nignar  ordar  flltar  tnat  utilizes 
tha  aacond  proeadura  by  cascading  2  or  aora 
sections  of  the  flltar  network  shown  In  Fig.  6.a. 
Tna  result  Is  a  nigh  ordar  fully  prograaaabla 
general  purpose  flltar,  that  can  ba  tailored  to 
aateh  alaost  any  propoaad  apectflestlon. 
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COMPUTER  SINULATIOHS  AND  EXPERIMENTAL 
VERlFlCATIOaSi 


Fig.  T  anowa  dlffsranot  eoaputar  slaulatlon 
outputs  of  tna  prograaaabla  flltar.  Tna  plots 
slaulata  tna  filter  raapoaaaa  asaualng  Ideal  OAa 
with  inflnlta  Cain  Banawtdgh  Products  (CBWP),  as 
wall  as  practical  'itar  raapenaaa  aaauaing  OA's 
finite  CBUPs  of  <  aa  gf  that  of  the  LM7AI  OA. 
A  single  pole  Oi  icial  aga  utilized  to  approzl- 
aata  tna  flltar  tr  itsrgr  flBaatlana  In  tna  later 
oasa.  The  sppr  uz.qMIm  was  found  adequate 
sines  tna  slaulu  >  jn  MMsitS  of  the  non  ideal 
rasponsa  ware  fou.id  ••  be  af  close  prozlalty  to 
the  siparlaantal  rsaulla  at  rig.  B.  The  aaparl- 
aantal  rasulta  ware  aatalaad  using  a  three  bit 
word  for  filter  topalagy  prograaaabl  1 1  ty  to 
aalaet  the  typo  of  tranafar  functions.  A  two 
words,  four  aits  aacn,  ware  used  for  filter  para- 
actars  prograaaabl  1  Ity  wnara  ■>.  and  Q.  ara 
controlled  Indapandantly  as  given  in  Tatla  3. 
Fig.  7  also  llluatrataa  a  nignar  ordar  prograa* 
aabttlty  wnara  a  fourth  order  cnaractaristlcs 
ara  shown  for  a  LPF  and  a  Chabycnev  BPF. 

S.  CONCLUSION 

The  novel  design  daaerlbad  hare  has 
resulted  In  a  universal  prograaaaola  filter  tnat 
can  ba  digitally  controlled  to  realize  alaost 
any  practical  flltar  spaclfleatlona.  This  la 
dona  through  the  use  of  CMOS  awltchaa  controlled 
by  binary  codas  to  prograa  the  ordar  of  tna 
filter,  the  filter  topology,  tha  filter  canter 
frequency  and  selectivity.  The  design  proeadura 
required  developing  optlaua  switching  arrange- 
asnts  for  tha  alnlnus  redundancy  in  components 
and  tha  least  dapandanee  of  tha  filtering 
function  on  swliening  laporreetlons  such  aa 
switches  stray  capacitances  and  non-zero  and 
non-ltnear  swltch-on  resistance.  Further  Invest¬ 
igation  Is  being  eonduetad  to  develop  a  program- 
oable  switched  capacitor  realization  that  can 
allow  frequency  scaling  by  changing  clock 
frequency.  Work  Is  also  In  progress  for 
developing  an  sztendsd  bandwidth  progrannabls 
filter  using  the  oomposlte  operstlonal  amplifier 
technique  proposed  earlier  by  the  author.  Such 
Implaosntatlon  would  lead  to  a  vary  useful 
aonollthio  devloe  at  moderate  cost. 
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